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Abstract 



This paper presents an assessment of global economic energy potentials for all major natural energy resources. This work is 
' ' based on both an extensive literature review and calculations using natural resource assessment data. Economic potentials are 
presented in the form of cost-supply curves, in terms of energy flows for renewable energy sources, or fixed amounts for fossil 
and nuclear resources, with strong emphasis on uncertainty, using a consistent methodology that allow direct comparisons to be 
made. In order to interpolate through available resource assessment data and associated uncertainty, a theoretical framework and a 
computational methodology are given based on statistical properties of different types of resources, justified empirically by the data, 
and used throughout. This work aims to provide a global database for natural energy resources ready to integrate into models of 
C*~) energy systems, enabling to introduce at the same time uncertainty over natural resource assessments. The supplementary material 
provides theoretical details and tables of data and parameters that enable this extensive database to be adapted to a variety of energy 
. . systems modelling frameworks. 
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1. Introduction 

Energy policy decisions for the planning of new energy gen- 
eration capacity designed to respond to future demand require 
information regarding the engineering feasibility and the re- 
quirements of such systems in terms of capital investments, but 
also, crucially, the availability of natural resources. Meanwhile, 
future energy systems are expected, in many contemporary pol- 
icy frameworks, to evolve towards their gradual decarbonisa- 
tion, in order to decrease anthropogenic interference with the 
climate system (see for instance Edenhofer et al. |U]). From an 
energy perspective, the decarbonisation of the sector involves a 
transfer of supply from traditional fossil fuel based technologies 
towards low GHG emission energy generation capacity such as 
renewable energy systems or nuclear reactors. Such a transfer 
requires changes in the technologies used through substitution 
processes, a subject extensively studied in the past (for instance 
I12j-|4j]). However, these transformations also require changes in 
the use of primary energy sources. Realistic energy scenarios 
of the future can only be designed in a way that does not ex- 
ceed natural sources and flows of energy which are available in 
all regions of the world. Therefore, assessments of the potential 
of natural energy resources are essential to energy planning and 
policy. 

Meanwhile, many models exist that generate scenarios for 
the evolution of global energy systems in the future (For a re- 
cent brief review of IMAGE/TIMER, MERGE, E3MG, POLES 
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and REMIND see Edenhofer et al. |Q]]). Such models, in or- 
der to generate scenarios that are realisable, must take into ac- 
count the limits to each type of natural resource. However, 
while some energy models do not currently take explicit ac- 
count of limits to resource flows, most of them do not con- 
sider their associated uncertainty. Of particular interest in this 
work is the Energy-Environment-Economy Model at the Global 
level (E3MG), a macroeconometric model of the global econ- 
omy, which calculates economic activity in 42 industrial sec- 
tors within 20 regions of the world . In a previous paper, a 
new sub-component designed for use in E3MG was introduced, 
modelling the trend of investor behaviour in the power sector 
facing the choice of various technologies, using their Levelised 
Cost Of Electricity production (LCOE), which includes the cost 
of natural resources, and the resulting technological substitu- 
tions In this model, the limitation of natural resources 
occurs through the use of cost-supply curves, requiring a com- 
plete set for each type of natural resource considered in order 
to properly constrain the model. 

Many studies and reviews have been published that sum- 
marise what was known of global energy potentials at time of 
their publication These reports however do not pro- 

vide any economic structure to energy potentials. Without un- 
derlying individual cost structures associated with energy re- 
source extraction or clarifications the use of the concept of 
economic potentials, natural resource limit values are ambigu- 
ous. As the key for strategic energy planning lies precisely with 
the cost structure of energy production, such assessments are 
of limited use for energy systems modelling and policy-making 
without additional assumptions over their cost distributions. 

Every energy source is limited, either in its total amount that 



can be consumed, for stock resources, or in the total energy 
flow it can produce at any one time, for renewable resources. 
Resources tend to be exploited in order of their cost of extrac- 
tion. As consumption gradually progresses to higher and higher 
levels of exploitation of particular resources, additional units 
consumed tend to incur increases in production cost. There- 
fore, the economic potential is better defined as a function of 
cost, rather than as a constant value. As the costs of production 
increase with the levels of use, developers increasingly seek al- 
ternatives, where they exist, through an evaluation of the op- 
portunity cost. Thus, the economic potentials for all energy 
resources available within a particular market depend onto one 
another and cannot be determined individually without know- 
ing the alternatives. Additionally, since every particular market 
for energy is composed in a particular way, economic potentials 
vary geographically. As resource use evolves and depletion pro- 
gresses, all economic potential values in a market change with 
time, an effect that stems from changing individual costs of en- 
ergy production, and and that has repercussions onto the rate of 
consumption of every natural resource. 

As discussed in previous work, the economics of energy re- 
sources when used, for instance, for electricity production can 
be modelled in a simplified manner using a complete ensem- 
ble of cost-supply curves, which express the cost of resources 
at various levels of exploitation [|4J]. In such a framework, the 
marginal cost of electricity production for every individual nat- 
ural energy resource using specific power technologies may be 
compared using a framework such as the LCOE at every level 
of natural resource use, in order to compare their profitabil ity. 
In such a model, the consumption of particular resources is lim- 
ited by increasing marginal costs, and the potential depletion of 
resources becomes naturally represented since the cost gradu- 
ally diverges at a levels of exploitation closer and closer to the 
technical potential. Cost-supply curves for global resources of 
wind, solar and biomass energy have been calculated previously 
using the land use model IMAGE 2.2 I18H22I1 . and are used in 
this way in the TIMER energy sub-model. Additionally, global 
cost curves for fossil fuels have been produced by Rogner [23 ], 
an influential work which is unfortunately becoming increas- 
ingly outdated. However, no comprehensive global assessment 
of all major energy resources which provides an underlying cost 
structure currently exists in the literature. 

Additionally, assessments of natural energy resources inher- 
ently possess high uncertainties, which must be taken into ac- 
count in order to generate confidence levels in model outputs. 
For instance, the global bioenergy potential has been estimated 
to lie between around 310 to 660 EJ/y, by Hoogwijk et al. i20ll . 
between and 650 EJ/y by Wolf et al. Q and between 370 
and 1550 EJ/y by Smeets et al. B25I1 . These particular uncer- 
tainties stem from those on future projections of food demand 
and levels of technology advancement in the agricultural sec- 
tor used in these calculations. These types of uncertainties are 
present in all assessments of renewable energy resources. Sim- 
ilarly, uncertainty arises with knowledge on amounts of stock 
resources, such as uranium and fossil fuels, where lower and 
lower levels of confidence are associated with larger and larger 
quantities. These uncertainties originate directly from the cu- 



mulative amount of effort that has been deployed to explore ge- 
ological occurrences, and express the fact that it is never possi- 
ble to know with certainty the detailed composition of the crust 
of the Earth. Using a review of literature, it is possible, and ap- 
propriate, to define ranges of energy potentials rather than spe- 
cific values, therefore defining areas in the cost-quantity plane 
for where actual cost-supply curves are likely to lie. 

This work proposes a theoretical framework and a compu- 
tational methodology for building natural resource assessments 
readily useable in models of energy systems, by using a com- 
bination of cost-supply curves and a treatment of uncertainty. 
This methodology is then applied to produce a cost-quantity 
analysis for every major natural energy resource, those with a 
potential larger than 10 EJ/y. As part of this work, cost-supply 
curves were produced for 13 types of resources for every one 
of the 20 world regions specified in E3MG, and form a new 
sub-model for natural resource use and depletion. Underly- 
ing potentials were however defined for 190 individual coun- 
tries, and can be aggregated to any other particular set of world 
regions. For the sake of presentation in this paper, the cost- 
supply curves were aggregated into global curves, providing a 
world view. Since data for all 190 countries could not realis- 
tically be provided here, tables of parameters are provided in 
the supplementary material that enable to reconstruct the cost- 
supply curves for a set of 14 world regions that were assumed 
the most representative of the requirements of the international 
modelling community. Additional resource specific informa- 
tion regarding theoretical derivations, additional methodology 
and justifications are also provided. 

This work follows consistently a methodology that can be 
reused as presented by the modelling community. To this end, 
definitions of the concepts used are first given followed by a 
concise description of the approach, detailing the cost-quantity 
analysis of resources with the associated uncertainty. Following 
this, a theoretical characterisation of the statistical properties of 
resource occurrences is given in order to enable the use of func- 
tional forms for interpolating resource data. This methodology 
is then used to produce cost-supply curves for renewables and 
stock resources. A summary of all major energy resources is 
given in the last section. 

2. Definitions 

Renewable and stock energy resources: Natural sources of 
energy that may be found in one of two forms: stocks, 
where energy may be extracted from fixed amounts of ge- 
ologically occurring materials with specific calorific con- 
tents; renewable flows, where energy may be extracted 
from continuously producing onshore or offshore surface 
areas with wind, solar irradiation, plant growth, river 
flows, waves, tides or various forms of heat flows. 

Theoretical Potential: Total quantity of energy stock or flow 
estimated to exist or stem from a particular natural process , 
disregarding its technical recoverability. 

Technical potential: Total quantity of energy stock or flow 
estimated to exist or stem from a particular natural pro- 
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cess, recoverable using a specific technique, disregarding 
its level of technical difficulty and the associated costs. 

Economic potential: Quantity of energy stock or flow esti- 
mated to exist or stem from a particular natural process, 
recoverable at exploitation costs that are competitive com- 
pared to all other alternative ways of producing the same 
energy carrier. Since the cost value considered competi- 
tive at any one time changes continuously, the economic 
potential is expressed as a quantity of stock or flow of en- 
ergy function of cost. However, this concept is used more 
conveniently in models when expressed as the inverse of 
this function, the cost-supply curve. 

Cost-supply curve: Function of the cost of energy flow or 
stock from a particular resource given that a certain quan- 
tity is already in exploitation or has already been con- 
sumed (the marginal cost). In this work, the cost-supply 
curve and the economic potential are used interchange- 
ably. 

Uncertainty range: Area in the cost-quantity or cost-flow 
planes in which actual real cost-supply curves have a 96% 
probability of lying. This would correspond to two stan- 
dard deviations, 2<x, if the distributions were normal, but 
they are in general skewed. Real values have a 2% proba- 
bility of occurring below the range, and a 2% probability 
of lying above. 1 

Productivity: Amount of energy stock or flow produced per 
unit of land or sea surface area (wind, solar, biomass, wave 
energy), bore depth or digging effort (oil, gas, coal, nu- 
clear fuels, geothermal energy) or construction effort (hy- 
droelectricity, tidal energy). 

Hierarchical resource distribution: Statistical type of natu- 
ral resource productivity distribution, with productivity 
values that strongly depend on the number of simultaneous 
positively contributing physical factors. Resource produc- 
ing items of this statistical type within one kind of resource 
(windy areas, rivers, mines, wells, etc) possess widely dif- 
ferent productivities which enable their ranking in order of 
resource quality. 

Distribution for nearly identical resources: Statistical type 
of natural resource productivity distribution in which re- 
source producing items possess nearly identical produc- 
tivity values, which do not depend on the simultaneous 
occurrence of several factors. Producing items of this type 
within one kind of resource (for instance sunny or fertile 
plots of land) have nearly identical properties, cannot be 
ranked and can be exchanged for one another. 

3. Methodology 

3.1. Econotnic potentials 

Natural occurrences of energy resources are found in differ- 
ent forms, with varying productivity levels or require various 



levels of effort for their extraction, which enable their transfor- 
mation into usable energy carriers with different levels of prof- 
itability. These variations together lead to particular distribu- 
tions of costs for their utilisation. Naturally, resources with the 
best qualities for energy production, and thus lowest extraction 
costs, are likely to be considered first by energy firms under 
financial constraints. Therefore, deriving economic potentials 
for energy resources involves the task of classifying and rank- 
ing different occurrences of specific resources in order of cost. 

Information on energy resources is scarce and irregularly dis- 
tributed, possibly inconsistent, and thus must be organised and 
classified in order to produce a consistent and complete set 
of economic potentials. Data may be patchy and incomplete, 
in which case assumptions are required in order to interpolate 
through missing parts. Such assumptions are taken in this work 
in the form of functional forms for the ranking of resources in 
terms of their cost of extraction. These are derived theoretically 
from basic statistical properties of resources. They have been 
carefully verified against several sets of data for specific types 
of natural resources which do not take any assumptions over 
the distribution of resources (wind, solar, two types of biomass 
resources as well as with uranium). They have been assumed to 
hold true for all other types of resources (for fossil, geothermal, 
hydroelectric and ocean resources). 

Methods of assessment differ and produce different results 
or ranges of results. In the absence of justified criteria onto 
which to base a choice of particular studies over all others, re- 
source assessments must be considered equally, the collection 
of which can be used to generate uncertainty ranges. This al- 
lows to decouple this work from specific assumptions used in 
specific studies. 

The methodology presented here builds upon the approach 
defined in our earlier work 0. The economic potential of re- 
sources is defined using the cost-supply curve, which expresses 
the quantity of resource available for any cost value consid- 
ered economic, or competitive with all other alternatives. Such 
curves are derived from cost rankings of resources and resulting 
distributions. The cost variable, however, stems from varying 
levels of technical difficulty for extracting resources, or alterna- 
tively, the productivity of energy producing resources such as 
plots of land, mines, oil wells, rivers, etc. Therefore, continuous 
distribution functions for the amounts of resources available in 
nature are defined in terms of their productivity. Two empirical 
forms for these distributions are defined and used throughout. 
Confidence ranges are derived from uncertainty analysis. The 
combination of both is used to construct probability densities 
for the location in the cost-quantity or cost-flow planes where 
the real cost-supply curves would be situated if it were possible 
to determine them with certainty. 2 These probability distribu- 
tions may be used as inputs to uncertainty analysis (such as 
Monte Carlo simulations) in energy systems modelling. 

3.2. Cost-supply curves 

The ranking of resources in terms of their productivity, re- 
quired for building cost-supply curves, can be done using a set 
of histograms of the quantity of energy stocks of flows that can 
be obtained within various ranges of productivity values. The 
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Figure 1 : a) Sketch of a hypothetical distribution of cost ranked amounts of energy or energy flow units available in various cost ranges. The uncertainty over the 
amount available in each cost range is indicated with a colour shading: the top of the dark histogram represents the minimum amount which has a 98% probability 
of being exceeded while the top of the white histogram indicates the maximum amount associated with a 2% probability of being exceeded. The most likely 
quantity is intermediate, represented by the top of the grey histogram, b) Cumulative distribution of energy resources, with uncertainty shown as vertical error bars, 
c) The marginal cost, or cost of extracting an additional unit of energy stock or flow given that a certain quantity has already been exploited, commonly called the 
cost-supply curve, with uncertainty represented as horizontal error bars, d) Cost-supply curve defined as a probability distribution, where the red curves indicate 
the limits of a 96% confidence level region in the cost-quantity plane, while the blue curve corresponds to the most probable cost-supply curve. The assumption is 
therefore taken that there is a 2% chance that the cost-supply curve lies below the upper boundary, and a 2% chance that it lies below the lower boundary. 



productivity variables may be converted into costs, which re- 
sult in a new set of histograms representing the amount of en- 
ergy that can be produced at costs within various cost ranges. 
This is shown in Figure Q] a), with a typical distribution of en- 
ergy resources, which decreases at low cost values due to the 
decreasing number of resources of exceptional quality, and to 
high cost values due to the decreasing productivity or recov- 
erability of the resources. The shading is a representation of 
the confidence associated to their potential availability. The top 
of the dark grey distribution shows the lowest quantity of as- 
sured resources, assumed to be exceeded with a probability of 
98%. The top of the white histogram represents the upper range 
of speculative resources, those assumed to be exceeded with a 
probability of 2%. The quantity which is the most likely to be 
available lies between these extremes, shown with the top of the 
grey histogram. 

The amounts in each cost range thus possess an uncertainty. 3 
In order to determine the quantity that can be obtained at or be- 
low certain cost values, the cumulative distribution function is 
calculated, shown in panel b). This sum converges towards the 
technical potential. The uncertainty increases approximatively 
cumulatively with increasing cost values through the root of the 
cumulative sum of the squares of the individual uncertainty val- 
ues, shown with blue error bars. 

The marginal cost of resources, or the cost of extracting an 



additional unit of resource given that a certain quantity has al- 
ready been used, corresponds to the inverse of the cumulative 
sum, shown in panel c). Thus, the cost of additional units di- 
verges when the number of units used approaches the technical 
potential, at the point of resource depletion. Using the uncer- 
tainty ranges, or error bars, to define two additional curves, as- 
sumed to define the upper and lower 96% confidence limits, a 
probability density can be defined in cost-quantity space for the 
location where the real cost-supply curve would lie if it were 
possible to know it with certainty. This is shown in panel d), 
where the red curves define the uncertainty area and the blue 
curve represents the most probable of all possible cost-supply 
curves, the mode of the distribution. Such probability densities 
are normally skewed, since the uncertainty over undiscovered 
resources lies at higher quantity values. 4 Note that the uncer- 
tainty is assumed to vanish at the contemporary position in the 
cost-quantity plane, since current costs and levels of exploita- 
tion are well known. 

3.3. Distributions 

Natural resources are scattered around the planet in different 
forms with different probabilities for the cost of their exploita- 
tion. Complex processes underlie the formation of these dis- 
tributions, however, they may follow certain statistical trends, 
the nature of which stems from the nature of the resource. One 
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Figure 2: Depiction of two types of natural resources, based on their statistical 
properties. Top. A typical sharp distribution for nearly identical resource types 
is shown in green, while the broad blue distribution is for hierarchical resources, 
from equations ^ and (2)- Both are expressed in terms of productivity. Differ- 
ent colour shadings represent uncertainty, as in panel a) of Figure \l\ Bottom. 
Same distributions expressed as functions of cost, from equations j5J and |4) 
through equation l|3j. Associated cost-supply curves are given in the inset. Note 
that the technical potential was adjusted to be the same for both curves for vi- 
sual clarity. 



particular property affects significantly these statistical distri- 
butions, whether the productivity of resource producing units 
(rivers, plots of land, wells, mines etc) tend to be very sim- 
ilar, which makes their ranking difficult, as opposed to re- 
source types which have strong ordering. For example, solar 
energy can be produced using photovoltaic (PV) panels, and 
these panels may be installed with equal ease almost anywhere, 
and scarcely populated regions of similar solar irradiation will 
have large potentials for solar energy situated within very nar- 
row productivity ranges where every plot of land has nearly 
the same productivity value. These are described by a particu- 
lar type of statistical distribution for nearly identical resources. 
The properties of wind energy potentials follow a very different 
structure. Plots of land within a geographical region are un- 
likely to possess the same average wind speed, and can there- 
fore be ranked in order of productivity, and thus of cost for the 



production of wind power. As opposed to solar energy sites, 
wind power production sites can be described by a hierarchical 
statistical distribution type. 

Resource distributions for nearly identical resources are 
sharply defined in a narrow range of productivity, but cut off at 
a maximum value, which corresponds to the best possible con- 
ditions for energy production. Meanwhile, resources with hier- 
archical distributions occur in large numbers in low productiv- 
ity ranges, and in ever lower amounts as productivity increases. 
This property stems from the large number of positively con- 
tributing factors which are required simultaneously in order to 
have a resource producing unit with high productivity. Such 
a property results in a distribution that decreases exponentially 
with increasing productivity, but cut off below a certain low 
productivity value, where it is simply assumed that no energy 
can be obtained with any reasonable amount of effort. This is 
shown in Figure [2] top panel, where a typical distribution for 
nearly identical resources is shown in green, and a hierarchical 
distribution is given in blue. Similarly to panel a) of Figure Q] 
the colour shading indicates uncertainty. 

These resource distributions are well described by the fol- 
lowing density functions: 
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where v is the productivity (see definition above), A is the tech- 
nical potential, cr is the width of the distribution and p,\ is the 
minimum usable productivity, in the first case, and fi2 is the 
maximum productivity available in the second case. 

Costs are related to the productivity by an inverse relation- 
ship, 

C = ^ + C , (3) 
v 

with which the distributions can be transformed into cost- 
quantity space. The scaling factor C var corresponds to a cost 
per unit of land or sea area (wind, solar, biomass, wave energy), 
digging effort (oil, gas, coal, nuclear fuels, geothermal energy) 
or construction effort (hydroelectricity, tidal energy), and the 
ratio C var lv has units of cost per unit of energy produced. Co 
corresponds to the sum of costs which do not depend on the pro- 
ductivity, per unit of energy produced. The conversion of these 
productivity distributions into cost distributions is described in 
detail in section S.2 of the supplementary material, and yields 
the following: 
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where A is the technical potential, B a scaling factor and Co a 
cost offset, the set of three parameters required to define every 
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distribution given in this work. These functions are illustrated 
in Figure [2] where the inset shows the associated cost-supply 
curves. It is observed that for a similar technical potential, the 
curve for nearly identical resources possesses less curvature up 
to very near the technical potential than those for hierarchical 
resources, a property that stems from their lack of ordering , and 
results in similar cost values for most of the resources. 

These functional forms have been found to reproduce very 
closely the cost-supply curves calculated by Hoogwijk et al. 
using the land use model IMAGE, whose work does not as- 
sume any functional dependence on cost for its distributions 
lfl8 . T9L 221. Distributions were calculated by producing cost 
ranked histograms of calculated potential renewable energy 
flows (wind, solar and biomass) at every point of a 0.5°x0.5° 
grid of the earth's onshore land. Thus, their form originates 
purely from statistical properties of the aggregation and rank- 
ing of the resources modelled. Using least-squares non-linear 
fits, the cost-supply curves in their work were found to agree 
very well with one or the other of the functions given above, 
depending on the nature of the resources: solar energy and agri- 
cultural land are well represented by the distribution for nearly 
identical resources, while wind power and rest land are well 
represented by the hierarchical distribution. Additionally, the 
distribution for hierarchical resources was found to agree well 
with observed cost distributions of uranium as reported by the 
International Atomic Energy Agency (IAEA) II26I1 . Examples 
of non-linear fits of these functions to IMAGE data are given in 
section S.2.5 of the supplementary material. 

No such global cost ranked data exist for the remaining types 
of natural resources that could enable fits of distributions dis- 
tribution. Choices of distributions were therefore taken as as- 
sumptions, made based on the physical nature of the resources. 
Potential basins that could be created for hydroelectricity pos- 
sess very individual characteristics, which makes them hierar- 
chical. Geothermal resources, however, were treated as a hy- 
brid mixture of the two, since good geothermal sources in ac- 
tive volcanic areas such as Iceland can be ranked, but large 
amounts of very similar sites can be found in non active ar- 
eas. Stock resources however were treated slightly differently. 
Different oil and gas occurrence subtypes (e.g. conventional 
gas, shale gas, clathrates, etc) originate from different geolog- 
ical processes which have no strong relation to one another, 
and should therefore be treated independently. These resource 
subtypes are characterised with different costs of extraction. 
They were assigned one hierarchical distribution per subtype, 
and subsequently aggregated. This resulted in composite cost- 
supply curves with complex structures. Coal, uranium and tho- 
rium resources were considered to occur as a single subtype, 
supported by the fact that the data for uranium was found to 
follow well the hierarchical distribution. 

3.4. Uncertainty 

The methodology used in this work for treating uncertainty is 
fundamental to this analysis of economic potentials as it allows 
the incorporation all available information, even when sources 
are inconsistent or conflicting. Inconsistencies can be found 
between assessments for most individual natural resources, and 



are the result of the use of different approaches and assump- 
tions, which can be determinant for the technical potential val- 
ues derived. This is most obvious in resources such as wind 
power, solar energy and bioenergy, where the total amount 
of appropriate land depends highly onto competing activities, 
making the assumptions in the evaluation of the land suitability 
factors the main drivers of uncertainty. Other such assumptions 
are world population and the associated food demand, levels of 
technological development and changes in agricultural produc- 
tivity. Resource assessments are uncertain by nature, since it 
is not possible to know with certainty the complete geological 
content of the crust of the earth, or to predict the weather and 
associated wind, sunshine and rainfall with perfect foresight. 
Thus the comparison of ever larger numbers of natural resource 
assessments is the key to define ranges of confidence, and these 
are as important as their associated most probable potential. 

This work uses a consistent methodology to define proba- 
bility distributions for cost-supply curves. Three cost-supply 
curves are derived from resource assessment data, where two 
are used to define the 96% confidence region in the cost- 
quantity plane, and one taken as the most probable of all pos- 
sible curves. In all plots of this work, the most probable cost- 
supply curves are given in blue and the 96% confidence lim- 
its are displayed in red. Uncertainty ranges are almost always 
asymmetric, since upper ranges are intrinsically characterised 
by smaller amounts of accumulated knowledge. 

Uncertainty is treated differently for renewable resources 
compared to stock resources. For renewables, cost-supply 
curves were obtained from the literature or calculated and taken 
as the most probable curves, while the 96% confidence limits 
were obtained by scaling the technical potential to the limits of 
its uncertainty range, defined by an ensemble of different stud- 
ies. 5 In the case of stock resources, all resource assessments 
provide classifications associated with cost ranges and various 
levels of confidence. In these cases, three cost-supply curves 
were calculated by assigning probabilities to uncertainty classi- 
fications, as in panel a) of FigureQ](i.e. reserves were assumed 
to exist with a 98% probability, while reserves plus all specula- 
tive resources were assumed to be available with a 2% proba- 
bility). Individual methodologies for all types of resources are 
described in the supplementary material. 

4. Renewable energy resources 

4. 1 . Wind energy 

Wind speeds depend strongly on altitude as well as on land- 
scape topologies, the climate and the type of land cover, or 
roughness. In general, wind speeds increase logarithmically 
with elevation at low altitude (see for instance S0rensen I127ln . 
and, for a specific elevation and geographic location, occur sta- 
tistically following a well defined Weibull probability distribu- 
tion which decreases both towards zero and large wind speeds 
(for instance Grubb and Meyer ll28lo . Average wind speeds on 
sites useable for energy production, for instance in the United 
Kingdom, range between 5.1 and 9 m/s at 10m elevation, the 
lower boundary determined by technology and the upper limit 
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Figure 3: Cost-supply curves for renewable resources: wind power, solar energy, hydroelectricity, geothermal power, biomass and ocean energy. 
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by the decreasing supply of locations with large wind speeds 
(for instance in Il27l0 . The power density offered by land ar- 
eas must be calculated using technical characteristics of partic- 
ular turbines. For a particular site, one integrates over all wind 
speeds the product of the site wind speed probability distribu- 
tion and the turbine power curve. However, a correlation exists 
between the yearly averaged wind speed and the number of full 
load hours [29]. The minimum possible distance between a tur- 
bine and its neighbours in a wind farm are determined by losses 
produced by the wake of neighbouring turbines which results 
in lower wind speeds and increased turbulence, and scales with 
the turbine rotor diameter, limiting the density of energy that 
can be extracted per unit of land area (for instance in the work 
of Mackay 1 30]). 6 Hoogwijk et al. assumed a maximum density 
of energy production of 4 MW/km z 1 19]. 

Various research groups have calculated the global distribu- 
tion of wind power, resulting in a range of values between 70 
and 2509 EJ/y for onshore wind power lll9[ |28, 31-34|, and 
about 57.4 EJ/y for offshore wind power [16]. A global on- 
shore value of 346 EJ/y has been derived Hoogwijk et al. fl^l . 
in whose work, used for the present analysis, estimations of av- 
erage wind speeds were applied to points on a global onshore 
grid, as well as the land suitability for the installation of wind 
farms using the land use model IMAGE 2. 2. 7 Energy potentials 
obtained from yearly averaged wind speeds determined on ev- 
ery point of the grid were subsequently aggregated into various 
cost-supply curves for specific regions of the world, according 
to the land suitability factor of each point. Cost values were de- 
termined using a present value calculation including fixed and 
variable costs, capacity factors and energy densities associated 
to particular land areas. 

Wind farm sites follow very strongly the distribution for hi- 
erarchical resources (see for instance the exceptional fit of Fig- 
ure S.4.1 in the supplementary material), since good sites with 
average wind speeds exceptionally suitable for energy produc- 
tion are geographically scattered, and the majority of areas in 
any region of the world possess mediocre average wind speeds, 
allowing a strict ordering of resource units (see for instance the 
European Wind Atlas lf35ll ). The profitability of a wind farm 
venture depends strongly on the quality of the wind resource, 
determined through capacity factor and average turbulence val- 
ues. For a fixed turbine investment cost, low capacity factors 
increase dramatically the cost per unit of electricity produced. 

Figure [3] presents the global economic potential for wind 
power. It gives an aggregate cost-supply curve, using for the 
most probable curve the data calculated by Hoogwijk et al., 
which involves the most detailed methodology for determin- 
ing the suitability factors, resulting in a technical potential of 
346 EJ/y 1 19]. While Lu et al. estimated an optimistic techni- 
cal potential of 2509 EJ/y by calculating wind potentials over 
the global onshore area and excluding low wind areas by re- 
stricting capacity factors to values above 20% [31], Archer et 
al. calculated a potential of 2257 EJ/y in an assessment where 
the land included was restricted to class 3 wind energy sites 
but did not include alternate uses of the land, a value taken as 
the upper boundary of the uncertainty range ll33ll . Meanwhile, 
an estimate of 70 EJ/y was obtained using an evaluation of the 



number of sites with average wind speeds above 5.1 m/s, but 
with an arbitrarily chosen value of 4% of that land available 
for wind turbine installation, in order to account for alternative 
land uses, taken as the lower boundary of the uncertainty range 
1 34 ]. All other existing studies result in values within this range 

liim. 8 

It is to be noted that using up a large fraction of that po- 
tential results in large areas becoming covered by wind farms. 
Since typical individual wind turbines currently have capacities 
of 3 MW but capacity factors of about 25-35%, compared to 
80% for coal power stations, replacing one coal power station 
of 1 GW for wind energy requires 700 to 1 100 turbines, cover- 
ing an area of 500 to 1500 km 2 , compared to about 1 km 2 for the 
original power station 9 . Even though agricultural land used by 
wind farms may still be cultivated, and therefore a competition 
for land with agriculture does not directly occur, strong emis- 
sions reductions pathways based on substituting fossil fuels for 
renewables likely implies that, given the large number of wind 
turbines required, these would invade permanently traditional 
rural landscapes. 

In the case of offshore wind power, although distributions of 
wind speeds at offshore locations tend to have a higher median, 
the air flow possesses less turbulence and wind speeds vary less 
in time, all of which contribute to produce a higher power den- 
sity in terms of geographical area, the total area where such tur- 
bines can be installed is small compared onshore areas, unless 
floating turbines become widely available [36]. The potential 
of offshore wind energy was evaluated to 22 EJ/y by Hoogwijk 
et al. and to 57.4 by Krewitt et al., both based on the work of 
Fellows et al. with power density values of 10 and 12 MW/km 2 
1 3 32, 37]. This potential could be at best only approximately 
six times lower than the most probable potential of onshore en- 
ergy given in Figure [3] and costs per unit energy are signif- 
icantly higher [38]. Due to the lack of reliable and consistent 
cost data, a cost-supply curve for offshore wind is not presented 
in this paper. 

4.2. Solar energy 

Solar radiation over the Earth surface is of about 
1.2xl0 5 TW, or 3.6xl0 6 EJ/y Q. The fraction of that en- 
ergy that can be harvested with existing systems has been es- 
timated by several studies, and ranges between about 1340 to 



14800 EJ/y |U8l|21l|40j]. Even though the total generation of en- 
ergy from solar technologies has been increasing steadily dur- 
ing the last two decades, they still represent very low percent- 
ages in their respective categories; solar heating systems ac- 
count for 0.3% of the total energy used for heating in 2008, and 
solar electricity generation represented only 0.06% of the total 
electricity generation during the same year [41]. 

Solar energy can be harvested using either of two existing 
technologies, photovoltaic (PV) devices (see B42I1 ) or concen- 
trated solar power (CSP) (see the International Energy Agency 
(IEA) Energy Technology Systems Analysis Programme (ET- 
SAP) [43]). Single crystal silicon photovoltaic diodes currently 
have light conversion efficiencies of up to 25%, while III-V 
semiconductor cells such as GaAs systems have efficiencies of 
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up to 28%, and solar cells using concentrated sunlight can con- 
vert light as efficiently as 43.5% [44]. The resulting electricity 
generation energy density ranges between 5 and 100 MW/km 2 , 
depending on the type of devices used and the geographical lo- 
cation. On the other hand, CSP technology uses a traditional 
steam turbine where water was heated using sunlight concen- 
trated with various arrangements of mirrors, and have efficien- 
cies of around 13-24% and energy densities near 25 MW/km 2 
HUH]. These systems are however restricted to high irradi- 
ance areas, and therefore have a lower global technical poten- 
tial. Solar energy is well represented by a distribution for nearly 
identical resources, since within areas of similar irradiance and 
average cloud coverage, its cost does not depend strongly on 
the nature of the land, having identical productivity values that 
depend solely on the chosen technology, and the opportunity 
cost of the land is the limiting factor to its technical potential. 

Figure [3] shows the global economic potential for solar en- 
ergy, using PV as technology. It is an aggregation of curves de- 
termined in various world regions, based on the work of Hoog- 
wijk 111 811 . but rescaled to match the values of de Vries et al. 
[21]. The global technical potential is of 3384 EJ/y. Regional 
cost-supply curves were drawn from an analysis performed us- 
ing the land use model IMAGE 2.2 to determine land suitability 
factors and the opportunity cost of the land at every point of a 
global grid, while regional estimates of solar irradiation were 
used to determine the energy potential, which have been taken 
here for the most probable cost-supply curve. The lower bound- 
ary curve of the uncertainty range, with a technical potential of 
1340 EJ/y, is the one evaluated by Hoogwijk et al., while the 
most probable curve was taken as the 2000 value of de Vries et 
al.. The upper boundary curve of the uncertainty range, with a 
technical potential of 14778 EJ/y, is the 2050 projection of de 
Vries et al., and stems mostly from increases in land availabil- 
ity following future reductions in the amounts of land required 
for agriculture through improvements in productivity I12ll . 

The use of a large fraction of the technical potential for solar 
energy using PV systems signifies covering up large amounts 
of land with solar panels. Mackay calculates an average pro- 
ductivity value of 10 MW/km 2 for the United Kingdom, while 
Hoogwijk used values between 6 and 25 MW/km 2 , depending 
on the geographical location, with capacity factors of up to 50% 



the extent that around 23% of the global hydroelectric technical 
potential is currently in use. However, its exploitation around 
the world is not even: 25% of the European technical potential 
has already been developed, while Africa uses only 7% of its 
hydroelectric resource l!46ll . 



[18 



30]]. When compared to coal power plants of capacity of 
1 GW and capacity factor of 80%, this implies that the replace- 
ment of one such power plant by solar panels requires an area 
between 50 and 500 km 2 . 

4.3. Hydroelectricity 

Hydropower stems from water pressure gradients that are 
produced by the run-off of rainfall through landscape topogra- 
phies, using dams to restrict water flow and accumulate water at 
elevation level higher than that given by the landscape, produc- 
ing a potential for electricity production using turbines (see for 
instance J45ll ). Hydroelectricity is the most deployed renewable 
electricity technology, with a global installed capacity of close 
to 1 TW, which produced around 2% of the total primary energy 
supply in 2008 114 111 . As a fraction of its total technical poten- 
tial, it is also the most developed of all renewable resources, to 



Costs of hydroelectric systems are highly site-specific 
and were found to have varied between around 400 to 
4500 USD2002/kW in an extensive global analysis done by 
Lako et al. M47I1 . These values are influenced by many differ- 
ent factors, which include material and labour costs, but also 
critically the opportunity cost of the land. The latter refers 
to the consequences of flooding large areas of land, and the 
resulting displacement of communities and agricultural activ- 
ities, and thus varies strongly from region to region. For this 
reason, the deployment of hydropower is often decided on po- 
litical rather than financial grounds. Hydroelectric resources 
were assumed in this work to follow a hierarchical distribution, 
since available natural basins that can be flooded possess vastly 
different geographic characteristics that make them unique and 
produces strong ordering. 

Figure [3] shows the global cost-supply curve for hydroelec- 
tricity, where the dotted line indicates the cost at the current 
deployment of 12 EJ/y, 23% of the modest value of the most 
probable technical potential of 66 EJ/y, calculated in this work 
from data gathered in the World Atlas and Industry Guide 201 1 
i46ll . The high deployment to potential ratio is an indication 
that the remaining number of suitable sites for building dams 
is relatively limited. The intersection of this curve with the 
current total hydroelectricity generation value yields a cost of 
production of about 68 USD/MWh. However, the development 
of hydropower projects hardly follows an order based onto cost, 
but follows instead an order dictated by political considerations, 
which are out of the scope of this work. Therefore, this value is 
only indicative, and projects with LCOE values between 23 and 
460 USD/MWh have been recently developed H. The use of 
this cost-supply curve in modelling involves the inevitable as- 
sumption of development following a cost order. In long term 
scenarios, this is reasonable, since the development of the lim- 
ited number of remaining available sites involves either increas- 
ing opportunity costs in inhabited areas due to increasing lo- 
cal populations, or increasingly large transmission costs associ- 
ated with increasing distances to uninhabited areas. The cost- 
supply curve was derived using the theoretical, technical and 
economic local potential values from j46ll . Since the definition 
of the economic potential in IJHD is not given, the (asymmet- 
ric) range of the distribution of recent cost values in the work 
of Lako et al. was interpreted (in 2008 dollars) as what is cur- 
rently assumed economic. The remaining technical potential 
(above the economic potential) was assumed to involve higher 
costs. The upper boundary curve of the uncertainty range was 
derived by considering the aggregated global theoretical poten- 
tial of 144 EJ/y from the data in the World Atlas and Industry 
Guide, while the lower boundary curve of the uncertainty range 
was derived by assuming that no additional construction of hy- 
droelectric dams occurs in the future, limiting future hydroelec- 
tric generation to 12 EJ/y. 
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4.4. Geothermal energy 

Geothermal resources, stored beneath the Earth's surface in 
the form of heat, are heat sources constantly replenished by 
the radioactive decay of isotopes of uranium, potassium and 
thorium (see for instance H49l l50l0 . Although geothermal heat 



has been used since prehistory, and its utilization for electricity 
generation commenced at the beginning of the last century, its 
current deployment is small in comparison with other sources 
of energy. It currently accounts only for 0.3% and 4% of the 
total electricity generation and heating production respectively 

Geothermal resources are classified in four categories: hy- 
drothermal (liquid and vapour dominated), hot dry rock (where 
fluids are not produced spontaneously), magma (molten rock in 
regions of recent volcanic activity) and geopressured (hot high- 
pressure brines containing dissolved methane) ll53ll . The most 
commonly used type is hydrothermal, although high expecta- 
tions exist regarding the development of Enhanced Geothermal 
Systems (EGS), oriented towards the hot dry rock type through 
hydraulic stimulation 15411 . According to estimations made by 
Aldrich et al. based on a report of the Electric Power Research 
Institute (EPRI), the estimated geothermal heat accumulated 
under the continental masses to a depth of 5 km depth is of 
approximately 1.46xl0 8 EJ, most of it assumed to be stored in 
rocks and water, with a proportion of 6: 1 in favour of the former 
1 55 , 56| . Even though this is a vast amount of heat, only a small 
part is recoverable for productive purposes. 10 

While geothermal resources are available all over the world, 
their accessibility differs from site to site according to vari- 
ous technical characteristics including the geological structure 
of the ground and the depth and type of heat reservoirs. In 
the vicinity of tectonic plate boundaries, narrow zones char- 
acterised by significant volcanic activity (so-called volcanic 
belts), geothermal gradients are particularly high, between 40 
and 80 °C per km of depth, enabling the extraction of high 
temperature geothermal resources. On the other hand, areas 
with low volcanic activity are characterized by low and uni- 
form geothermal gradients: around 25 °C per km of depth 
lf55ll5al . The extraction of geothermal resources in active areas 
are highly site-specific, and thus were assumed to follow a dis- 
tribution for hierarchical resources [p7l . Meanwhile, geother- 
mal gradients in the rest of land masses have very similar prop- 
erties and costs, and were therefore assumed to follow a distri- 
bution for nearly identical resources. Cost-supply curves were 
produced for both types of land and both hydrothermal and EGS 
technologies, generating four curves which were subsequently 
aggregated in each world region. 

Stefansson found a high correlation between the number of 
active volcanoes in a particular region and the estimate of the 
size of hydrothermal resources for electricity generation in the 
same region 115 811 . Using the total number of volcanoes active in 
the world, discarding those located on the sea floor or in arctic 
regions, he estimated a global installable hydrothermal electric- 
ity producing capacity of approximately 200 GW (producing 
6.0 EJ/y of electricity with 95% capacity factor). Using this in- 
formation, along with the statistical analysis between wet and 



dry systems developed by Goldstein et al. 115911 . Bertani esti- 
mated the total geothermal installable electricity producing ca- 
pacity of 1200 GW (36 EJ), including hydrothermal and EGS 
technologies 

The global aggregation of curves yields the cost-supply curve 
presented in Figure[3] The associated global technical potential 
of 36 EJ/yr, involves a 95% capacity factor. Cost values were 
obtained from the International Energy Agency (IE A) [62]. The 
lower and upper boundaries of the uncertainty range, of 4 and 
114 EJ/y, are explained in section S.3.4 of the the supplemen- 
tary material. 

4.5. Bioenergy 

Bioenergy, energy derived from plants, is currently the most 
widely exploited renewable energy resource, with 51 EJ/y, 10% 
of global annual primary energy use 114 ill . The combustion 
of biomass derived fuels is nearly carbon neutral if CO2 up- 
take during plant growth is taken into account, minus losses 
occurring in transformation processes. Thus, biomass based 
technologies provide an important emissions mitigation poten- 
tial. While biomass combustion using integrated gasification 
combined cycle technology (BIGCC) is expected to become 
the most efficient biomass based electricity production method 
1 63], the combination of biomass and carbon capture and stor- 
age technology has been shown to produce negative CO2 emis- 
sions 1 64], thus providing the potential for reductions of atmo- 
spheric CO2 concentrations, or for compensating other emis- 
sions. Moreover, the emissions factors of some power plants 
using conventional coal technologies are being be reduced by 
co-firing coal and biomass fuels. Meanwhile, liquid biofuels 
derived from biomass, such as ethanol and biodiesel, have the 
potential to replace oil-derived transport fuels with minimal 
changes in vehicle internal combustion engine technology and 
jet engines l65ll . 

Biomass currently used for electricity and biofuel produc- 
tion largely originates from forestry and agricultural residues, 
and other forms of commercial or household mixed solid waste. 
Volumes of waste available could amount up to 100 EJ/y but 
are highly uncertain and not studied here 11 I25ll66ll . The larger 
share of bioenergy potential lies with the production of dedi- 
cated biomass crops. Global technical potentials for primary 
bioenergy range between and 1550 EJ/y l20l 12 ll l24l l25l l67ll . 
Bioenergy crops include perennial woody short rotation cop- 
piced trees, such as willow, poplar or eucalyptus, perennial 
grasses such as miscanthus, elephant grass and switchgrass, 
starch rich crops such as wheat, corn, sugar beet and cane, and 
oil rich crops such as rapeseed and palm. Depending on their 
nature, they can be transformed into energy carriers by using, 
among many processes, combustion, gasification or anaerobic 
decay for electricity production, fermentation or the Fischer- 
Tropsch process for transport fuel production ll68l . 

Biomass production for energy purposes makes use of agri- 
cultural land and thus may have a high opportunity cost. The 
technical potential that lies in agricultural land is large, but en- 
ergy production from biomass is in direct competition for land 
with food production, a situation which has the potential to 
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drive significant increases in world food prices [69]. Follow- 
ing the approach which has been used recently by many au- 
thors (for instance in |20[ 22, 24, 25]), the explicit assumption 
is taken in the present work that future bioenergy production 
uses no more than leftover land after the global food demand 
has been met, a premise that is difficult to justify in the absence 
of specific legislation and further investigations, but it avoids 
the complex problem of simulating food and biomass prices. 12 
Thus the bioenergy potential is obtained by subtracting from 
the total biomass potential the amount required by the food de- 
mand, based on population growth curves and dietary assump- 
tions. 

Hoogwijk et al. evaluated the use of the land at each point 
of a global grid yearly up to year 2100 using IMAGE 2.2, in 
which leftover agricultural land was termed 'abandoned land' 
1I20L 221. The reported cost-supply curves were observed in 



the present work to follow a distribution for nearly identical 
resource units using non-linear fits of eq. (|5} to their data. In 
addition to agricultural land, however, other types of geograph- 
ical areas with lower productivities exist which can be used for 
particular bioenergy crops. These were labelled 'rest land' by 
Hoogwijk et al. and contribute a significant global technical po- 
tential. They were found to follow the distribution for hierar- 
chical resources by using fits of eq. © to their data. Examples 
of such fits are described in section S.2.5 of the supplementary 
material. Land use depends strongly on assumptions regarding 
world population, diet habits, global urbanisation and trade of 
agricultural products. The four main SRES scenarios, Al, A2, 
Bl and B2 (see ll70Tl ). were taken as assumptions for all exoge- 
nous variables in these calculations, and results are presented 
for each. Large differences arise between scenarios, with tech- 
nical potentials ranging between around 302 EJ/y for the A2 
scenario to 676 EJ/y for the Al scenario, which result in large 
uncertainties for values of the global biomass technical poten- 
tial. This work has however been revisited more recently taking 
into account additional assumptions of future water scarcity and 
land degradation by van Vuuren et al., yielding lower estimates 
ranging from 120 to over 300 EJ [67]. Other ranges have been 
estimated using different methodologies, with more pessimistic 
projections of to 648 EJ/y by Wolf et al., and optimistic values 
of 367 to 1548 EJ/y by Smeets et al. JUS]. The low end of 
the range given by Wolf et al. stems from high projected food 
demand and low agricultural productivity, while the high end 
is due to mostly vegetarian diets and high productivity. Mean- 
while, the high end of the range of Smeets et al. originates from 
'super high' agricultural productivity, high availability of the 
land and landless animal production systems. 

Figure [3]presents the global economic potential of bioenergy 
in terms of primary energy before conversion to electricity or 
liquid biofuels (derived from the data in using 
both abandoned and rest land, but includes also a small com- 
ponent from bagasse of 3 EJ/y. Four cost-supply curves are 
given, calculated by Hoogwijk et al. for the Al, A2, Bl and 
B2 SRES scenarios in 2050, shown as solid curves 13 M22I1 . A 
value near zero was taken for the lower boundary of the un- 
certainty range, consistent with the low end of the range calcu- 
lated by Wolf et al. while the high end of the range calculated 



by Smeets et al. was taken for the upper boundary of the un- 
certainty range [24, 2511 . For a decarbonisation scenario, the 
cost-supply curve derived for the Bl SRES scenario was con- 
sidered the most probable cost-supply curve, but for other types 
of scenarios, choices of curves consistent with particular work- 
ing assumptions should be made. 

4.6. Ocean energy sources 

The term ocean energy denotes renewable energy produced 
using seawater as a resource, where unlike for wind energy or 
hydroelectricity, not only the kinetic energy of seawater can be 
used to produce electricity, but also temperature gradients in the 
ocean and salinity differences near river mouths. Using ocean 
energy as a general classification type, it can be divided into 
four main sources of energy 117 U 17211 : 



• Wave energy, driven by transfers of energy from the wind 
to the surface of the ocean, 

• Tidal energy, driven by the rise and fall of sea levels due 
to gravitational forces (tidal range) and the resulting water 
currents, 

• Ocean Thermal Energy, driven by temperature gradients 
between upper and lower ocean layers, 

• Salinity Gradient energy, derived from salinity gradients 
between ocean and fresh water at the mouths of rivers. 

Section S.3.6 of the supplementary material provides a re- 
view of theoretical potentials for these sources, resulting in a 
total that could be as high as 600 EJ/y. Technical potentials 
however are much lower and uncertain, since the current devel- 
opment status for ocean energy technologies excluding tidal is 
preliminary, and cost data is in some cases unavailable. Spe- 
cific geographical and configurational requirements for tidal 
and salinity gradient technologies involves, as it is the case for 
hydroelectricity, calculating the technical potential by summing 
the potential values from a large number of individual studies. 
Such studies have not been performed exhaustively on a wide 
scale yet. Meanwhile, wave and ocean thermal are based onto 
global extrapolations carried out using physical measurements. 
Global energy potentials calculated in various studies are given 
in table [1] and additional details are given section S.3.6 of the 
supplementary material. 

Energy potentials for ocean thermal and salinity gradient en- 
ergy are theoretical and highly uncertain, and no reliable cost 
estimates were found. These types of resources were therefore 
not included in the present calculations for the most probable 
cost-supply curve, due to the risk of generating misleadingly 
optimistic potentials given the lack of reliable information. 14 
Wave and tidal systems are better established. Therefore, us- 
ing cost values obtained from the IEA and the ETSAR a cost- 
supply curve for ocean energy based on an aggregation of sepa- 
rate cost-supply curves for wave and tidal energy was produced, 
and is given in Figure [3] with a small technical potential of 
22.5 EJ/y [48, 71]. The lower and upper boundaries of the un- 
certainty range were obtained from the extremal values of 8 and 
72 EJ/yr respectively given in table [1] 
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Figure 4: Cost-supply curves for fossil and nuclear resources, including oil, gas, hard coal, soft coal, uranium and thorium. Hard coal includes anthracite and 
bituminous coal, defined as coal with a calorific content above 16 500 kJ/kg. Soft coal corresponds to sub-bituminous coal and lignite, and includes all coal with a 
calorific content lower than 16 500 kJ/kg. 
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Table 1 : Technical potentials for different types of ocean energy used to define 
the cost-supply curve. The uncertainty ranges are defined using the Min and 
Max values, while Mode represents the most probable value. 



5. Stock resources 

5.7. Fossil fuels 

As it occurs with all types of exhaustible natural resources, 
fossil fuel resources and reserves are known to continuously ex- 
pand, even though they are gradually consumed. This is due to 
periodic resource discoveries and improvements in the methods 
of extraction. Therefore, what is considered economical to ex- 
tract changes every year. Reserves are distinct from resources, 
the former referring to the resources that are known to exist 
with almost complete certainty and to be economical to extract, 
while the latter refers to those which are thought to exist with 
various degrees of confidence, and those currently thought too 
expensive to extract. As technological improvements and addi- 
tional knowledge affect the economics of different methods of 
extraction, there is a flow from resources towards reserves, and 
thus reserves expand fl23l l79ll . Meanwhile, discoveries continu- 
ously add to resources. As prospecting activities for hydrocar- 
bon sources remains very active, this makes the production of 
cost-supply curves more difficult than for renewables, since is 
at best a snapshot in time of what is known to exist and recov- 
erable with current technology. 

In order to assess global energy potentials, it is nevertheless 
necessary to explore cost-supply curves for fossil fuels, even 
if they are derived from current knowledge, and therefore ex- 
pected to change in the future. It is unlikely that fossil fuel 
resources turn out smaller than what is currently expected to ex- 
ist. On the contrary, it is probable that they turn out significantly 
larger as methods of extraction are devised for types of occur- 
rences which were until recently not thought possible to use, 
such as gas hydrates or oil shales. Therefore, the cost-supply 
curve uncertainty ranges are highly asymmetric. The associ- 
ated extraction costs, which increase as low cost conventional 
sources are depleted, nevertheless decrease due to technologi- 
cal improvements, and it is therefore not immediately obvious 
whether costs are likely to go up or down in the future. 

Global cost-supply curves have been calculated previously 
by Rogner (23J. These results have been used extensively by 



the energy modelling community; however they are becom- 
ing increasingly outdated. This section provides an update to 
the work of Rogner, but using an approach emphasising uncer- 
tainty, and thus, following the spirit of the current treatment of 
renewable resources, in opposition to the approach of Rogner, 
the results of this section should be interpreted as ranges rather 
than specific values. 

The economic potentials of fossil fuels are given in Figure |4) 
showing in order liquid hydrocarbons, gaseous hydrocarbons, 
hard and soft coal, the last two being classified using their 
calorific content. 15 For oil and gas, different types of occur- 
rences considered in this assessment are indicated with text. 
These are associated with independent distributions of the hier- 
archical type, aggregated to produce composite curves. Due to 
the wide use and global diffusion of fossil fuel extraction tech- 
nology, extraction cost ranges were assumed to be the same for 
all regions of the world. In the case of coal, less information 
was found over differing types of mines and associated costs, 
and single distributions were used, where costs were assumed 
to vary little with the amount extracted. This is unlikely to mat- 
ter in the long run given the very large scale of the resource, and 
limited expectations of its depletion. 

Uncertainty ranges were determined using resource classifi- 
cations, and in some cases where this is unavailable, their na- 
ture. Oil occurrences obtained from the World Energy Council 
and the German Federal Institute for Geosciences and Natural 
Resources are classified as either reserves or resources, with 
the exception of oil shales, which are as a whole considered 
resources only IU3U15I1 . Four types of oil resources were con- 
sidered, conventional (crude) oil, oil shales, oil sands and extra- 
heavy oil. Cost ranges were obtained from the IEA [80]. Gas 
occurrences follow a similar trend, but a larger number of types 
of resources were considered: conventional gas [13], shale gas 
18111 . ti ght gas 111 ill 1311 . coalbed methane B82ll and methane hy- 
drates fl83[] . The respective cost ranges were obtained from the 
ETSAP 18411 . Large amounts of methane are known to exist dis- 
solved in aquifers ll ill , but were not included due to the lack of 
reliable data. Information for coal was derived from a mixture 



of data B13l 11511 . Complete details on the methodology under- 
lying these curves, as well as region specific data tables can be 
found in the supplementary material. 

5.2. Fissile materials 

Five sources of fissile materials for nuclear reactors are 
known to exist. These are enumerated in order of cost. The 
first comes from stocks of highly concentrated 235 U (uranium) 
or 239 Pu (plutonium) originating from decommissioned nuclear 
arsenals diluted with 238 U. The second source is lightly en- 
riched 235 U/ 238 U produced from mined natural deposits. The 
third originates from U and Pu recovered from spent fuel (using 
the PUREX process B85ll ). The fourth source is thorium (Th) 
using the 232 Th/ 233 U fuel cycle. The fifth source is U which oc- 
curs in very low concentrations in seawater. Producing a cost- 
supply curve involves creating a scenario for the nuclear sector, 
and requires careful consideration of uncertainty. Additionally, 
if ingenious use of fast reactors is invoked for the future, fuel 
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efficiencies of up to 50 times larger could be obtained, altering 
dramatically these expectations. 

In order to construct a cost-supply curve for U and Pu, the 
nuclear industry was assumed to continue to use current meth- 
ods and thermal reactors, and therefore, only fuels originating 
from naturally mined U and from nuclear arsenals were con- 
sidered. Many authors stress that deposits of Th worldwide are 
three times larger than those of U [85 88 1. However, less ef- 
forts at prospecting for Th ore have been carried out and as a 
consequence, the current reasonably assured reserves of Th, in 
tonnes of natural Th, are lower than those of U ll26ll . a situation 
which is likely to change if interest in Th grows. The nuclear 
fuel cycle for Th being more efficient than that of thermal reac- 
tors based on U, it leads to larger amounts of energy per tonne 
of natural Th and thus leads to lower fuel costs per unit of en- 
ergy. Costs only include the extraction costs given by the IAEA, 
without the inclusion of enrichment or transformation compo- 
nents. 

Detailed resource data from the IAEA for naturally occur- 
ring U and Th were used to construct two cost-supply curves 
and associated uncertainty ranges II26I1 . The data are classi- 
fied into four levels of certainty and four cost ranges. Such 
resources generally increase naturally in size with increasing 
costs of extraction, as well as with uncertainty, an effect pro- 
duced by the hierarchical ordering of natural resource consump- 
tion and by the decreasing amount of effort which has been 
spent on prospecting activities for resources more and more dif- 
ficult to exploit. For the conversion of resources from tonnages 
to energy values, an average conversion efficiency for thermal 
U reactors of 159 TJ/t was used, determined from the 2008 
electricity production of 2611.1 TWh from a global fleet ca- 
pacity of 273.7 GW, with a capacity factor of 80%, which used 
59 065 t of natural uranium [26]. Meanwhile, the burnup rate 
for Th reactors was derived from the value of 24 000 MWd/t 
reached by the experimental Indian model [86], equal to about 
2100 TJ/t. Panels e) and /) of Figure [4] present the resulting 
global economic potentials for U and Th. Uncertainty ranges 
for U were obtained by considering only reasonably assured re- 
serves (RAR) for the lower boundary of the uncertainty range, 
RAR and inferred reserves for the most probable cost-supply 
curve, and all of the RAR, inferred, prognosticated and specu- 
lative resources for the upper boundary of the uncertainty range. 

The uncertainty ranges are highly asymmetric due to the ten- 
dency of the size of speculative resources to increase with the 
level of uncertainty. It is observed that in terms of energy, re- 
serves of Th are larger than those of U, and that these are also 
less expensive per unit of energy, due to the higher burnup rate 
of the Th system. It must be emphasised that U resources could, 
in principle, be used with much higher bumup rates, were fast 
reactors to be deployed globally. The resources of U do not 
include seawater U, as data over these are scarce and highly 
speculative. 16 Finally, it is to be noted that the fuel component 
of the levelised cost for nuclear reactors is very small compared 
to the investment costs, which results in a very small influence 
of the fuel costs onto the decision-making, unless nuclear re- 
sources become depleted. 



6. Summary of energy resources 



Resource 


Use 


Technical Potential 


Name 


Type 


Dist. 


EJ/y 


L 


M 


U 


Units 


Wind 


Flow 


Hierarch. 


.72 


72 


346 


2257 


EJ/y 


Solar 


Flow 


Identical 


.04 


1340 


3384 


14778 


EJ/y 


Hydro 


Flow 


Hierarch. 


12 


12 


66 


148 


EJ/y 


Geotherm. 


Flow 


Hybrid 


0.23 


4 


36 


111 


EJ/y 


Biomass 


Flow 


Hybrid 


51 





447 


1548 


EJ/y 


Ocean 


Flow 


Hierarch. 


.002 


8 


23 


72 


EJ/y 


Oil 


Stock 


Hierarch. 


170 


9 


67 


98 


10 3 EJ 


Gas 


Stock 


Hierarch. 


109 


7 


46 


106 


10 3 EJ 


Hard Coal 


Stock 


Hierarch. 




24 


220 


419 


10 3 EJ 


Soft Coal 


Stock 


Hierarch. 


139 


5 


37 


75 


10 3 EJ 


Uranium 


Stock 


Hierarch. 


30 


0.83 


1.36 


3.43 


10 3 EJ 


Thorium 


Stock 


Hierarch. 




1.74 


4.68 


12.27 


10 3 EJ 



Table 2: Summary table for all energy resources. Stock/ Flow indicate whether 
resources are renewable flows or stocks. Hierarch. I Identical I Hybrid identi- 
fies the type of statistical distribution assigned. Use refers to current yearly 
consumption of these resources. L indicates the lower boundary of the uncer- 
tainty range. M indicates the most probable technical potential. U indicates the 
upper boundary of the uncertainty range. 



Table [2] provides a summary of all types of global energy 
resources, classified by type (renewable flows or stocks), to 
which a type of statistical distribution it is assigned (for hier- 
archical or nearly identical resources, or a hybrid mixture of 
both), along with technical potential values. The potential val- 
ues are given with their lower and upper boundaries of the un- 
certainty ranges. For comparison, current consumption of these 
resources is given based on data from the IEA MIA . Note that 
biomass resources are expressed in terms of primary energy, 
which become smaller when converted into electricity accord- 
ing to the efficiency of transformation. 

7. Conclusion 

This paper presents an assessment of global economic energy 
potentials for all major energy resources, those with a potential 
larger than 10 EJ/y. These were given in the form of cost-supply 
curves, adding an economic structure to energy potentials, and 
therefore providing them an unambiguous definition. Addition- 
ally, these were provided using a probabilistic construction that 
allows a simple representation of uncertainty. The curves were 
calculated using assumptions over the cost distribution of re- 
sources using functional forms based on statistical properties of 
resource types. The set of energy potentials include six types of 
renewable energy sources, wind, solar, hydroelectric, geo ther- 
mal, biomass and ocean energy, as well as four types of fossil 
fuels, oil, gas, hard and soft coal, and two nuclear materials, 
uranium and thorium. While the potentials for renewable re- 
sources were determined predominantly based onto an exten- 
sive review of the literature, potentials for stock resources were 
determined directly using resource and reserve assessment data. 

The cost-supply curves calculated in this work were pro- 
duced in order to be used by the global energy modelling com- 
munity, for the purpose of constraining models of the energy 
sector in order to produce realistic scenarios of future energy 
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use. It is hoped by the authors that this work will supersede out- 
dated studies currently used and provide a consistently calcu- 
lated update for all types of energy resources. In particular, the 
large set of regional cost-supply curves underlying the aggre- 
gate curves presented in this paper form the core of a new model 
for natural resource use and depletion for the global E3 model 
E3MG, to be used through the family of technology models 
FTT. Other regional aggregates can be provided by the authors. 

Resource assessments, however, change continuously as new 
information becomes available, and as new and more sophis- 
ticated studies are carried out. Therefore, while the potentials 
presented here may become outdated in the future, however the 
methodology presented will not. The simple and robust theoret- 
ical framework and computational methodology presented here 
can be very useful for reporting natural resource assessments, 
and enable direct use in models of energy systems. It would 
therefore be appropriate to reuse this approach with new data 
as they become available, either by rescaling the curves given 
(the values for the A parameters) or by recalculating new sets 
of parameters. In either cases, this approach provides a con- 
sistent and general framework for limiting all types of natural 
resources in energy models. 
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around 65 mW/m 2 at the continental level, producing an average thermal en- 
ergy recharge rate of about 315 EJ/yr [ 90] . This value can be considered as the 
theoretical potential of geothermal energy if viewed in terms of sustainable ex- 
traction of geofhermal resources over an extended period. However, the amount 
of time over which geothermal resources could be used at higher rates than this 
is likely to be more than one thousand years. 

1 1 Waste amounts depend on efficiency of biomass use (such as in food or 
timber production) and therefore subject to significant changes depending on 
future policy. These are therefore difficult to model, as are their associated 
costs of production. 

12 Bioenergy potentials could in principle be larger if global food demand 
is not met. However, it will not be lower if global food demand is indeed 
met. The problem of simulating food prices is complex as it involves modelling 
both local food markets, underreported in developing countries, and efficiency 
changes associated with changes in food prices. 

13 Values projected for the year 2050 were used since amounts of land poten- 
tially available in the future are expected to increase due to increased agricul- 
tural efficiency. 

14 The costs for ocean thermal and salinity gradient energy systems are likely 
to be much higher than those of tidal and wave energy. This would result in a 
piecewise cost-supply curve featuring an additional step at high costs. 

15 Hard coal includes anthracite and bituminous coal, while soft coal includes 
sub-bituminous coal and lignite, the last two having lower calorific contents 
than the first two. The limiting calorific value used to separate the two cate- 
gories is of 16 500 kj/kg Oil 

16 U is present dissolved at very low concentrations in seawater (3-4 ppb, 
from flS[ ). giving rise nevertheless to large amounts of U given the size of 
the terrestrial body of seawater. Water turnover due to currents is very slow 
however, making it highly speculative whether a significant portion of seawater 
U can be recovered, and the costs involved are very high [85]. 



Notes 

'The 2<x probability range correspond to erf( V2) = 95.45%, yielding 2.28% 
as a probability of values occurring above or below the range. The values of 
96% and 2% are used instead for convenience. 

2 Note that the use of uncertainty ranges in the cost-quantity/cost-flow plane 
relaxes the constraints of using specific functional forms, since it allows varia- 
tions in the particular forms of the functional dependences within the ranges. 

3 No error bars are present for the cost variables, since an uncertainty on 
costs corresponds to an uncertainty on how to distribute energy quantities be- 
tween existing cost ranges, which translates into an uncertainty in the quantity 
in each range. 

4 Thus, the most probable cost-supply curve is neither the mean or the me- 
dian of the skewed distribution, it is the mode, or maximum. 

5 This is done in order to avoid inconsistencies where curves calculated in- 
dependently, for instance by fitting data, could cross in some cases. 

6 While larger turbines intercept a larger wind front, they are also spaced 
further apart in two spatial directions. Thus, while the power production of 
large wind farms scales with the square of the length of the blades, it scales 
inversely with the square of the distance between turbines. These two effects 
almost cancel each other out, except for the fact that taller turbines intercept 
higher wind speeds at higher altitudes. 

7 For details on IMAGE see (§§]. 

8 Fellows et al. concludes with an estimate for 2020 of 148 EJ/y, while Grubb 
and Meyer calculated a global potential of 191 EJ/y I28j [32ll . 

'The variation originates from both assumed ranges in capacity factors of 
25-35% and turbine densities of 2.2 to 4 MW/km 2 GIUH]. 

10 Note that the average replenishment of the geothermal heat underground 
is several orders of magnitude inferior to the stock of heat currently available: 
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Supplementary Material 

S.l. Introduction and use of the supplementary material 

This part of the work aims to complement the main paper by providing all details that could be required by researchers interested 
in either: 



• Verifying the methodology or reproducing the results 

• Rebuilding this database for use in a particular model of energy systems 

A high amount of care was put into summarising compactly all relevant information in this part of the work in order to make this 
possible. Mathematical details underlying the calculations given in the main text are provided. Additionally, tables of data and 
parameters are given for a chosen set of world regions, which may not necessarily correspond to the particular divisions used by 
other research groups. It is however impractical to provide larger tables involving all countries of the world, even though such 
tables exist underlying this work (190 countries exist in our database). For more information, the authors may be contacted at the 
address provided. 

Natural resource assessments are performed continuously and what is known of global natural resources changes continuously. 
Therefore the cost-supply curves in this work may become outdated. This will happen due to three processes: firstly the total 
amount of resources might change (parameter A), secondly the scaling of costs may change through inflation (parameter B) and 
the costs of technology may reduce due to learning-by-doing (parameter Co). However, the structure of the cost-supply curves, or 
shape, will not change. These parameters should be simply scaled to new values and the results will still be valid. 



S.2. Derivation of distribution functions and cost-supply curves 

S.2.1. Distribution function for the hierarchical type of resources 

Hierarchical resources have an exponential energy distribution in productivity space: 

where A is the technical potential and cr is the half width of the function. This function is required in cost space, and the equation 
connecting cost C to productivity v is 

C = ^ + C„. (7) 
v 

where C var corresponds to costs per unit of effort or per unit of resource producing items such as the rent of the land (in $/km 2 ), 
bore hole depth ($/km), dam size and type, mine depth, etc. The ratio of C var to v has units of $/GJ. Co is the sum of fixed costs, 
such as capital investments, transport or transformation costs, etc, in $/GJ. 
The density productivity interval must be transformed into a cost interval: 

C C 

dC = r dv, ^dv = - dC. (8) 

v 1 (C - c y 

Using the cost scaling parameter B = C var /o~, the distribution becomes: 



AB 



■e c - c » dC C > C f ixed 



f(C)dC = \ (P-Cof -~ "'"J™ 8 . (9) 
C < Co 



S.2.2. Distribution function for nearly identical resources 

In the case of nearly identical resources, there is no convenient exact analytical form that can be derived from eq. |2]of the main 
paper. However, the form given in eq. |4]works extremely well when compared to data (see section [S.2.5b , and can be derived from 
eq.|2]through a simple approximation. 

Nearly identical energy producing resources, such as land plots, are assumed truly identical, and therefore have a potential 
situated at a single value of productivity v = //, 

n(v)dv = N6(v - fi)dv, (10) 
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where n is a density of energy producing land area, while N is the total energy producing land area (in km 2 ) and the function 6{v) 
is the Dirac delta function^ Without any additional reductions in productivity, the total amount of energy that can be obtained from 
these land resources would be their area times the productivity: 



A= f vNS(v — fx)dv = N/j.. (11) 
Jo 

Unit land areas have a suitability factor, however, that reduces their productivity below the maximum value of fi by a small amount 
e with a probability P. The probability for the reduction in productivity is assumed to be normally distributed around zero but 
positive e, with standard deviation much smaller than the average productivity <x << fi: 

P{e)de = { Vlb e ~^ e °^ e ^ , (12) 
\ otherwise 

where the reduction in productivity e must be less than the maximum value fi. The distribution of resources must be calculated by 
summing over all reduction values e given their probability P(e): 

f C" 2N _ « 2 
n'(v)dv— I n(v)P(e)dvde — I — e 6(v — n + e)dvde. (13) 

J J() V27TO- 

This is a convolution and can be seen as a sum of several Dirac Delta functions centred at slightly reduced values of productivity, 
fi - e, with probability P(e), instead of one Dirac Delta function centred at fi with probability 1 . The total amount of energy that 
can be obtained from each plot of land corresponds to its area times its productivity. Thus, the productivity distribution of energy 
production potential leads to eq. |2]of the main paper: 

< 2/v iv ~ t " 1 

- -ye ir? dv v < jj. 



g(v)dv=vn'(v)dv=i ^Tna . (14) 

{ v> n 

This function is required in cost space, and the equation connecting cost to productivity is 

C= — +C fixed , (15) 
v 

where C var corresponds to the rent of the land (in $/km 2 ), while Cfi xe d is the sum of fixed costs (in $/GJ). The productivity is situated 
very near the value of fi, since the variations of productivity are small and <x << jj.. v can be rewritten as a small variation around fi, 
i.e. [i - A: 

7 + L fixed - 1 T~, <~ Cfixed — H + Cfixed - — n 



C = + Cfi xe d = ^777 + Cf ixed =i — — | 1 + — | + Cf ixed = —J- (2yU - V) + Cfixed, (16) 



which is the crucial approximation, and using Co = C var /fi + Cfi xed , 

v = (Co (17) 

where Co is defined as the sum of fixed costs plus C var /fi, the total cost at the maximum productivity value. Since dC - -C var //idv, 
the density can be rewritten as 



gMcJ&fa-Q£ + *) eJS$tdc c>Co 

( C < Co 



(18) 

-o 

where the cost scaling parameter B is defined as C var o-//i z . This can be rewritten further as 

8 (C)dC = i Jb(^ + 1 K^ C C>C ° . (19) 
{ C < Co 



The Dirac delta function is defined such that S(x — a)f(x)dx = f(a). 
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The value of C cannot be below Co by definition, but also, the factor exp[-(C - Co) 2 /2£ 2 ] decreases rapidly to zero as C - Co 
becomes larger than B. Therefore, the value of the term (Co - C)/B is less than one wherever a significant potential of energy exists. 
Since cr is much smaller than p, this results with 

cr C - Co 



and therefore the distribution becomes 



« 1, (20) 



g(C)dC = \ l^B e 2B ~ dC C>C ° , (21) 
{ C < Co 

where now A = Np is approximately the total energy potential. 

This is the strict region of validity of the expression given in|4]of the paper. In numerical terms, it is empirically found that the 
rigidity of these rules can be relaxed and the validity extended. For example, the distribution in productivity space can actually have 
a tail towards higher values or have a large value for cr, and this does not significantly alter the goodness of fit of the function in 
cost space. 

S.2.3. Cost-supply curve expressions 

From the distributions f(C)dC and g(C)dC, cumulative distributions N(c) can be derived. For hierarchical resources, this results 

in 

JV(C) = Ai^), (22) 

while for nearly identical resources this is 

N(C) = A erf | (C ~ C()) |, (23) 
\ y/2B I 

where 'erf is the error function. 

The cost-supply curves are the inverse of these functions, which respectively give 

C(N) = + Co (24) 

and 

C(N) = V2B inverf (— ) + C , (25) 

where 'inverf is the inverse error function. 



S.2.4. Parameterisation formulas 

In order to define a particular distribution or cost-supply curve, the parameters A, B and Co are required, while the data available 
usually involves the total technical potential of the resource and a fraction of this considered to exist at costs situated between two 
values, as well as the current level of use of the resource. Assuming that two points of the curve are known, this may be expressed 
as two quantities Q\ and Qi which occur at two cost values C\ and C2. These quantities can be expressed as fractions of the total 
technical potential, 5\ and 62, the latter corresponding to the parameter A. In the case of the distribution for hierarchical resources, 
the values of B and Co are the following: 

Co = C2ln<?2 - Cllnt?1 , (26) 
In 62 - In 6\ 

B = -(Ci-C )\og6 1 . (27) 

In the case of nearly identical resources, this becomes 

fi=— , (28) 

V2(inverf<5i - inverf^) 

Co = V2B inverf 61 + Ci (29) 
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S.2.5. Demonstrating the validity of the functional forms using IMAGE data 

Examples of the use of the analytical forms of the distributions are presented in figure IS~4~Tl for biomass, solar and wind energy. 
The data are taken from land use simulations performed using IMAGE by Hoogwijk et al., which provide the only sources of 
cost-supply curves calculated outside of this project that do not already use assumptions on the analytical form of the resource 
distribution II 1 8ll 191 12211 . Since IMAGE simulates the use of the land on each point of a global grid, and since these cost-supply 
curves were calculated by building histograms of the number of grid points with productivities situated within various ranges, their 
form stems purely from the statistical nature of the data. These are thus appropriate for testing the functions given above. 

Non-linear least-squares fits were performed with both analytical forms for each data set. In every case, only one of the two 
functions given above is appropriate, while the other is not. Fits are moreover of exceptional quality. For instance, wind resources 
are the best example of resources of the hierarchically ordered type, which stems from the exponentially increasing number of 
simultaneous factors required to produce ever higher productivities. The data is found to follow almost exactly the distribution 
for hierarchical resources (note that the deviation at low cost values stems from the aggregation of a region with a different cost 
structure into the region for Canada). Meanwhile, solar resources represent the best example of nearly identical resources, since 
in regions of similar irradiation, all sun-facing areas are equivalent. The data is found to follow closely the distribution for nearly 
identical resources. Biomass resources from abandoned agricultural land are nearly identical. This stems from the similar nature 
of local areas of agricultural land (i.e. large plains, deltas, similar irradiation, etc). Land plots with lower productivity are used for 
other activities. Rest land, however, is the category of land which would not be used for agriculture, and can be of various nature, 
but includes mainly savannah, shrubland and grassland or steppe. These can be ordered, and can be seen to follow the distribution 
for hierarchical resources. 



S.3. Cost-supply curve calculation methodology per resource type 

S.3.1. Definition of world regions 

Cost-supply curves were calculated in this work for every E3MG world regions from aggregations of data defined for 190 
countries. However, the region definition in E3MG is very specific and does not correspond closely to that of most other global 
models, and tables provided here for E3MG regions would be of limited use to the global modelling community. For accuracy, 
data for 190 countries would be required to be provided here, but is not possible for space considerations. For the convenience of 
potential users, the results are provided in tables with a definition of regions resembling that of other models such as IMAGE, AIM, 
etc. Any other aggregation of data, in table or curve form, can be supplied by the authors upon request. Table IS. 4~T1 gives the list of 
regions used here with most countries that belong to them. 



S.3.2. Wind and solar energy 

Following the justification of section [S.2. 51 wind resources were modelled using a distribution of hierarchical type, while solar 
resources were modelled using a distribution for nearly identical resources. In both cases simulations performed by Hoogwijk et 
al. (wind energy), Hoogwijk (PhD thesis, wind, solar and biomass energy) and de Vries et al. (both) using IMAGE 2.2 were used, 
published in the form of data tables featuring both technical potentials and interpolations through cost supply curves at specific cost 
values for a list of 17 world regions II 1 8l 1 1 9l 12 1 11 . These values were used to find the distribution parameters A, B and Co for every 
one of their regions. In the case of wind, A values were thus obtained without additional processing. In the case of solar, A values 
were obtained from de Vries et al. while B and Co values from Hoogwijk IU8U21I1 . In both cases, B and Co values were obtained 
using equations l26lto[29l 

However, the regional aggregation in the aforementioned work does not match exactly the one chosen in the current study (or 
the one used in E3MG), detailed in section [S. 3. II In order to obtain curves for this set of world regions, energy potentials from 
IMAGE 2.2 regions were disaggregated into 190 countries, and subsequently re-aggregated. This required additional assumptions 
in particular cases where A values were required to be divided between underlying countriefl In the case of wind energy, the 
division of A values was done proportionally to the cube of the yearly averaged wind speecQ, times the amount of land suitable in 
each country for these energy production activitiesfthe land area times the suitability factor provided by Hoogwijk et al., assumed 
the same for all countries member of a region) !19lN93ll . In the case of solar energy, A values were divided proportionally to the 
insolation averaged over countries times the amount of land suitable in each country J91 , 93, 94 1. Assuming an identical shape for 
the cost supply curves (identical values of B and Co) for every country within a particular IMAGE region, and using A values thus 
divided, cost supply curves for the 190 countries were built. Given this set of curves, the re-aggregation of curves into new world 
regions was performed by summing the energy potential values at each cost (i.e. a sum along the horizontal axis of the cost-supply 
curve, called a horizontal sum henceforth). These aggregated curves do not correspond anymore to pure distributions of either type, 



3 Note that this is mostly true for E3MG regions; the regions used for this paper are very similar to those used by Hoogwijk et al. 1 19] 
4 Wind energy scales with the cube of the average wind speed averaged over time (see for instance [27]). 
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but do not differ significantly from pure forms in any of the regions chosen for this work. Thus, new values for A, B and Co for this 
work's regional definition were re-estimated using equations [26] to [29] for the sake of simple presentation in this work (avoiding 
listing parameters for 190 countries, or providing aggregate curves defined on large numbers of cost data points). For E3MG, data 
curves evaluated on 1000 cost data points are used directly instead^ 

Cost values with which cost-supply curves were calculated using equations l26l to l29l were also obtained from Hoogwijk et al., 



but were rescaled to 2008 prices H18L I19H . This procedure, however, generates costs of energy production slightly different (wind) 
or higher (solar) than recent estimates available from the International Energy Agency, due to small errors (wind) or significant 
learning-by-doing cost reductions (solar) stemming from economies of scale with large expansion of electricity generation capaci- 
ties that occurred between 2004 and 2008 [48]. The curves were therefore recalibrated with a constant offset to match recent values. 
The results are provided in table IS .4/3l for this work's list of world regions. 

S.3.3. Hydropower 

Hydroelectric resources, highly site dependent, were modelled using the distribution for hierarchical resources. Hydroelectric 
potentials and current annual electricity generation values were obtained from the last available technical report of the International 
Journal on Hydropower and Dams (IJHD), while the costs were obtained using an extensive study of 250 recent projects by Lako 
et al. from which statistics were derived [ 46l 47 1 . These statistics were performed for the countries studied in Lako et al., and were 



used as proxies for regions not studied in their work, or where no information on recent hydroelectric developments was found. 
Some countries do not have recently reported hydroelectric projects onto which to base cost values. 

Recent developments have hardly followed an order of cost, since they were scattered between 500 and 4000 2003USD/kW. 
In order to use a cost-supply curve, it can only be assumed that future developments actually will approximately follow a cost 
order. Although only approximately true, this is reasonable, since development costs will significantly increase when more and 
more usable sites are depleted, irrespective of the particular order in which they were built, and only difficult or distant river basins 
remain. This is important since the costs of hydroelectricity are currently not high in comparison to alternatives, but the resources 
are limited, and therefore the development of hydroelectric resources must be limited through an increasing cost in models of power 
systems such as FTT:Power. 

As can be seen in the current hydroelectricity generation data compared to the data for hydroelectric potentials from IJHD, a 
significant portion of the technical potential of every region is already developed lfl5l Pkill . The cost values delimiting the technical 
and economic potentials amongst remaining potential hydroelectric sites are not given by IJHD. Since the distribution of costs is not 
symmetrical, the assumption was taken that the amount of resources considered economical lies at costs between the local average 
cost ju minus its standard deviation 6, fx - 5, and plus twice its standard deviation fi + 26. This puts the upper cost limit to around 
5000 2008USD/kW. Thus, sites within the technical potential with costs higher than this are considered currently uneconomical. 
Given this definition, a cost-supply curve for each region was calculated. Parameters for each regional cost-supply curve are given 
in table IS. 4.31 The global cost-supply curve of figure [3] of the main text is an aggregation (a horizontal sum) of these regional 
curves. 

S.3.4. Geothermal energy 

Geothermal resources were divided into two groups, occurring in either "in belt" or "out of belt" land areas, referring to the so- 
called volcanic belt. "In belt" areas are located in volcanically active zones with high geothermal gradients (temperature gradients 
with bore depth from the surface of the ground). Given the particular characteristics of geothermal active areas in terms of their heat 
storage and underground temperature variation, the extraction of geothermal resources in those places are highly site-specific, and 
were therefore modelled using a hierarchical distribution. "Out of belt" areas corresponds to the rest of the continental masses, with 
sites that are characterised by smaller geothermal gradients, and that are almost identical to one another within large geographical 
areas. "Out of belt" resources were thus modelled using a distribution for nearly identical resources. The ratio of "in belt" to "out 
of belt" land area values were obtained from the 1978 report of the EPRI [56], enabling to divide reported technical potentials into 
two A values for each distribution type. Geothermal resources were moreover calculated for both hydrothermal and EGS dry rock 
technologies, yielding four sets of parameters. Each cost-supply curve in each region was obtained by aggregating four curves. 



Technical potentials for different world regions were obtained from Bertani 116111 . Given the differences between their regional 
aggregation and this work, the same methodology was used as for wind and solar energy in order to disaggregate the regional 
technical potentials between the same 190 countries. The proportion of the regional technical potentials assigned to every country 
within a particular region was assumed to be proportional to the total amount of geothermal energy stored up to five kilometres of 
depth in each country lf55ll . 

Cost values for geothermal electricity production were taken from the IEA [62]. It was assumed that 90% of the resources 
'in belt' were situated within these ranges of costs. However, resources 'out of belt' follow the distribution for nearly identical 



Exact analytical forms for cost-supply curves correspond to the inverse the cumulative distribution. When the cumulative distribution involves the sum of 
several distributions, an analytical form for the cost-supply curve does not exist. 
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resources, but face higher costs due to lower geothermal gradients and less productivity per unit investment. Since no additional 
cost information was available in this regard, these resources were assumed to lie in the upper half of the cost range given by [62]@ 
Table IS. 431 gives the parameters that can be used to reproduce these cost-supply curves using both types of distributions. 

The lower boundary curve of the uncertainty range assumes a technical potential of 4 EJ/y based on differing assumptions for 
both technologies. In the case of hydrothermal technology, a conservative potential estimate of 70 GW (2 EJ/y) was derived by 
limiting the calculation to well known sites that have been already characterised by direct involvement or informed calculations 
1 60, 61]. Meanwhile, the limited amount of accumulated experience with EGS technology creates uncertainties in the evaluation 
of the technical potential through variations in the efficiency of extraction 115411 . leading Bertani to estimate a lower limit of 70 GW 
(2 EJ/y) 16111 . The upper boundary of the uncertainty range involves yet another set of assumptions for both hydrothermal and EGS 
technologies. In the case of hydrothermal, according to Stefansson, undiscovered or additional resources could exist which would 
be five to ten times higher than identified resources, increasing the potential to 1000-2000 GW (57 EJ/y) lf58ll . In the case of EGS, 
the technical potential is calculated by an extrapolation of resources in the United States to the global level using the proportion 
between the heat stored at depths of less than 10 km in the United Stated with the known EGS primary energy potential in the same 
region, estimated as lxlO 6 EJ of heat stored per 2.61 EJ/yr of EGS primary energy potential [54]. Using the estimation of the heat 
stored at depths less than 10 km on the global scale of 403 xlO 6 EJ ll95ll . this estimation results in a global technical potential of 
57 EJ/yr, yielding a total of 1 14 EJ/y. 



S.3.5. Bioenergy 

Four bioenergy cost-supply curves are given in figure [3] of the main text, for each of the SRES scenarios Al, A2, Bl and 
B2, based on simulations performed using IMAGE 2.2 GamiH (see fzH for information on SRES scenarios). The primary 
biomass energy sources considered in the cost supply curve are abandoned agricultural land, rest land and bagasse, where the 
first is the largest source in all scenarios. Following the justification of section [S. 2. 51 abandoned agricultural land was modelled 
using distributions for nearly identical resources, while rest land was modelled using distributions for hierarchical resources. The 
remaining type of primary biomass resources, bagasse (from [15]), contributes very small fractions of the total potentials, and its 
technical potentials were simply added to the potentials from abandoned agricultural land and rest land, for every region in every 
scenario. Cost values, however, are only given by Hoogwijk et al. for the total amount of biomass resources in each region, not 
individually for abandoned agricultural and rest land [22] . Therefore, the right distribution to use had to be determined, by deciding 
which of the two represented best the data. This corresponds to finding the dominant distribution. Therefore, the appropriate type 
of distribution was determined by visual inspection for each region. Potentials for abandoned agricultural land are for most regions 
much larger than those for rest land, and therefore most regions were modelled using distributions for nearly identical resources. 
These distributions were disaggregated into 190 countries, following the methodology described in section IS. 3. 21 proportionally to 
country land areas times their suitability factor. Table I5~4~4| provides values that can be used to parameterise biomass cost-supply 
curves for the world regions used in this work, with the appropriate type of distribution used indicated in the last column. 



S.3.6. Ocean energy 

Given the vast extent of oceans, the calculation of theoretical potentials for ocean energy sources produces large values. For 
instance, using a global wind-wave model, Mork et al. estimated a potential for wave energy between 2986 and 3703 GWe (94 
to 117 EJ/yr), while Charlier and Justus estimated a global theoretical tidal power potential between 1000 and 3000 GWe (32 to 
95 EJ/yr) using a capacity factor of 100% IT751, l96ll . In the case of ocean thermal energy, Pelc and Fujita estimated a theoretical 
potential of approximately 10 TW (315 EJ/yr) using a capacity factor of 100%, while for salinity gradient energy, Cavanagh et al. 
calculated a value of 2.6 TW (82 EJ/yr) using a capacity factor of 100% J78, 97]. Using values from these particular studies, the 
total theoretical potential for ocean energy would be as high 523 and 609 EJ/yr. However, more reliable and conservative potentials 
have also been evaluated, given below. These values are modest in comparison. As indicated in the main text, cost-supply curves 
were calculated for wave and tidal systems only. For presentation in this work only, these two cost-supply curves were combined 
into a single one for ocean energy. Parameters for regional ocean cost-supply curves are given in table lS.4.3l 



Wave Energy 

In the case of wave energy, WEC estimated a maximum global installable capacity of 2 TW by limiting developments to tech- 
nically favourable locations near coastlines ll34ll . Using this value, and assuming a single capacity factor value of 32%, Rrewitt et 
al. estimated a technical potential for wave energy of 20 EJ/yr, while UNDP estimated a technical potential of 65 EJ/yr using the 
same value but assuming a capacity factor of 100% instead IH , T^. Following a more conservative approach restricted to shore- 
lines exceeding a power production of 30 kW/m (resulting in around 2% of global coastlines), Sims et al. estimated a technical 
potential of 500 GW (6.3 EJ/y using a capacity factor of 40%) B73I1 . It is clear however that using single capacity factor values is 



6 In 'out of belt' areas, the same technologies are involved, either for hydrothermal or EGS, as for 'in belt' areas. However, the resources are nearly identical over 
large areas and of equally low quality in comparison to 'in belt' areas. Significantly higher productivities are found in volcanic areas. 
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not appropriate. It is likely that, as it is the case for wind, the cost variation of a wave energy cost-supply curve should stem from 
capacity factor variations which stem from the local quality of the resource. Such data is however not currently available as it is for 
wind resources. Capacity factor distributions are likely to follow roughly those of wind energy, which vary between 15 and 35%, 
since both resources are closely related (wind capacity factor distributions were obtained by extracting the capacity factor from 
Hoogwijk et al. data fl9ll ). The assumption was therefore taken in this work that wave energy resources are captured using a single 
technology, with an investment cost given by ETSAP of 6600 USD/kW, with a capacity factor that varies between 35% (where the 
resource quality is highest, and the cost of electricity production is lowest per unit energy produced) to a low value of 15% (below 



which sites are not economically useable) 117 111 . Using a maximum global capacity of 2 TW, the cost-supply curves were calculated 
with a hierarchical distribution, assuming that 90% of the wave resources are available at capacity factors within the range 15-35%. 
The disaggregation into 190 countries was performed according to the lengths of their respective coastlines, using data from the 
Central Intelligence Agency (CIA) Q. 

Tidal Energy 

Accounting for most of the global installed capacity of ocean energy systems, tidal energy is the only technology that has reached 
a commercial scale, with approximately 523 MW installed at the end of 2010 Il98ll . WEC made a rough estimation of the technical 
potential of tidal energy of about 2000 TWh/y (7.2 EJ/y), 10% considered economical H S 1 1 0011 - In a more detailed study, 
Hammons presented a global but non-exhaustive list of potential tidal sites that could be considered for development, including 
projected installed capacities and approximate annual outputs B74I1 . The total output from these sites would be of almost 400 TWh/y 
(1.4 EJ/y). Hammons furthermore extrapolated that the inclusion of additional sites around the world not studied specifically in his 
work would result in a global technical potential for tidal energy likely to range between 500 and 1000 TWh/yr (1.8 - 3.6 EJ/y). 



The cost-supply curve for tidal energy was calculated using the range of cost values given in ETSAP of 5000 to 6500 USD/kW [71]. 
Existing capacity was assumed to have been built at costs below that range, while the sites reviewed by Hammons (400 TWh/y) 
were assumed to be associated with costs within the range [74]. Additional sites were assumed to have costs above the range. 

Ocean Thermal and Salinity 

The state of development of ocean thermal and salinity gradient energy technologies is currently experimental and therefore large 
uncertainties accompany calculations of associated energy potentials 07311 . Upper limits in the form of theoretical potentials have 
been calculated. Nihous estimated a theoretical potential for ocean thermal energy of 2.7 TW (85 EJ/yr or 23 652 TWh/yr), which 
corresponds to the maximum amount of energy resources that could be extracted without disrupting significantly the temperature 
of the upper layers of the ocean in an steady state regime, using a one-dimensional model of oceanic temperature gradients [76]. 
Using a similar method, Charlier and Justus produced a more conservative estimation for the theoretical potential of 1000 GWe 
(32 EJ/yr), assuming a capacity factor of 100% fl75ll . However, according to von Arx, such a level of heat extraction would imply a 



decrease in the ocean surface layer temperature of approximately 1°C II 1 111 . In order to avoid such a decrease, Charlier and Justus 
recommend a reduced estimate based on 10 TW of usable heat replenishment rate, corresponding to 100 GWe (3.2 EJ/y) IT75Tl . 

In the case of Salinity Gradient, based on average discharge and low flow discharge values, Skramesto estimated the theoretical 
potential in the range of 1600 - 1700 TWh/yr (5.8 - 6.1 EJ/yr) B77I1 . Using a global discharge rate of fresh water to seas of 44 500 km 3 
per year, Krewitt et al. estimated a theoretical potential of 2000 TWh (7.2 EJ/yr), value very similar to the estimate of Skramesto 



S.3.7. Oil 

Oil resources (Table IS.4.61 were considered in four types of occurrences, crude oil, oil shales, extra-heavy oil and oil sands 
|13lll5|]. Cost information was obtained from the IEA [80]. The data were aggregated into this work's world regions. For each type 



of occurrence for each region, a hierarchical distribution was parameterised by assuming that 1% of the resources have extraction 
cost below the lower bound, while 90% have a cost of extraction below the upper bound. The distributions were summed for 
each region in order to calculate regional cost-supply curves, and all distributions were summed in order to determine the global 
cost-supply curves given in figure |4]of the main text. 

The curve for the lower boundary of the uncertainty range was defined by assuming that only crude oil, extra-heavy oil and 
oil sands reserves are available, and that the rest is either unusable or does not exist. The most probable cost-supply curve was 
calculated assuming that crude oil, oil sands and extra-heavy reserves and resources are available, as well as oil shale resources, but 
no additional amounts. The curve for the upper boundary of the uncertainty range assumes that all reserves, resources and additional 
amounts are available, and that an additional amount of oil shales is discovered, evaluated at 50% of the current resources. This 
was done in order to compensate for the absence of speculative resources and lack of detailed information available for oil shale 
resources, which are likely to become larger if additional exploration is carried out, and will occur if (but only if) interest in oil 
shales intensifies 



7 For instance, if strong decarbonisation policies are implemented globally, oil shales are not likely to be explored much further, since they currently involve large 
processing costs and only small scale exploitation. 
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S.3.8. Natural gas 

Gas occurrences were considered in five forms, of which four unconventional: conventional gas, shale gas, tight gas, coalbed 
methane and methane hydrates fTH IT3I l8ll - [83ll . The associated cost ranges were obtained from the ETSAP [84]. Of the uncon- 
ventional forms, only shale gas has seen exploitation larger than experimental. An additional source of methane exists, which is 
thought very large, aquifer gas 111 ill . However, its potential being very speculative, no reliable information over costs of extraction 
was found, and thus these were not considered in the present study. Similarly, methane hydrates provide a very large source of natu- 
ral gas, however, these resources occurring under the sea, and the methods of extraction very experimental, the costs of exploitation 
are very large, and due to large amounts of shale gas available at lower costs, it is unclear whether the world will see wide-scale 
exploitation of methane hydrates. All resources were distributed into this work's world regions, except for the methane hydrates, for 
which it is not clear whether they are situated within territorial waters or not, and thus were assigned to an international category. 
Regional cost-supply curves were calculated with the same method as oil, and the global cost-supply curve is an aggregation of all 
regions. 

The curve for the lower boundary of the uncertainty range includes conventional gas reserves only. The most probable cost- 
supply curve includes conventional, shale and tight gas reserves, along with half the conventional, shale and tight gas resources. 
It moreover includes half of the coalbed methane reserves. The curve for the upper boundary of the uncertainty range includes all 
reserves and resources, including methane hydrates. 



S.3.9. Coal 

Although coal is a very common commodity and well known resource, information over its natural occurrences is not very 
detailed. Coal information was available from two sources, BGE and WEC, and table IS. 4. 81 was constructed using a mixture of 



both 113111511 . Where information was inconsistent, the larger amounts were kept (such inconsistencies were not frequent nor very 
large). Since BGR does not report coal resources in the complete classification (proven, probable and possible reserves or resources) 
for all countries, some elements of the table are nil 111 311 . This situation is likely to be due to the large amounts of coal available 
with conventional mining techniques, and therefore most of the resources are considered reserves, and occurrences with lower 
productivity or higher costs are not reported. Coal formations occur in different forms which have different calorific contents. For 
similar mining and transport costs, the costs of coal in terms of energy produced are higher for lower grade coals. Coal resources 
were divided into two categories, hard coal, including anthracite and bituminous coal which posseses higher calorific contents of 
between 16 500 and 35 000 kJ/t, and soft coal, including sub-bituminous coal and lignite, with calorific contents between about 
1 1 000 and 1 6 500 kJ/t lfl3ll . Note that these classifications are not strictly well defined in geological terms, and that coal occurrences 
exist that have intermediate properties. This stems from different geophysical processes taking place during the slow formation of 
these hydrocarbons. 

The curves for the lower boundary of the uncertainty ranges for soft and hard coal include proven reserves only. The most 
probable curves include proven and probable reserves, and half of the proven and probable resources. The curves for the upper 
boundary of the uncertainty ranges include all proven, probable and possible amounts for both reserves and resources. 



S.3.10. Uranium 

Information for uranium occurrences is available from a survey of the IAEA j26Tl . which compiles data provided by all member 
countries. They are highly detailed with four classifications and four cost ranges, as seen in table |S.4.9l As the amounts in each row 
of the table are cumulative (i.e. the data in one cell is inclusive of the sum of the cells to the left), with the associated cost values 
they correspond to cumulative distributions. Assuming that they should follow hierarchical distributions, the associated cumulative 
distribution may be fitted using a non-linear least squares method. Although a fit of a function with three parameters over four data 
points is hardly a reliable method to determine a best fit with any level of certainty, it nevertheless produced the best curves that 
could be interpolated between points, as was determined by close inspection of each fit. Note that additional data points were be 
defined in order to constrain the fits better, such as additional values at higher costs by repeating the last data point, assumed equal 
to the technical potential in the saturation region, and at (0,0). The resulting fits were found to follow the data very closely. Values 
for military stocks of U are uncertain, since the information that is publicly available is scarce and incomplete, and were omitted. 
These are very small in comparison to natural stocks [26]. 

The curve for the lower boundary of the uncertainty range includes only RAR (reasonably assured reserves) in all cost ranges. 
The most probable curve includes RAR and inferred reserves. The curve for the upper boundary of the uncertainty range includes 
all four classifications of resources in all cost ranges. The speculative resources in the unassigned cost range were not included in 
the non-linear fitting procedure, since their distribution into existing or higher cost ranges is ambiguous. They were therefore added 
to the technical potential of the upper boundary of the uncertainty range. Finally, it was assumed that sea water U is too costly and 
uncertain to include in the present study. Given the large amount of sea water on the planet, this resources is thought very large even 
though the concentration of U in sea water is very low. However, due to the very slow mixing process of sea water, it is misleading 
to consider the global body of sea water as a usable source of U, as it could be rapidly depleted locally, providing small amounts, 
without access to the remaining resources situated far offshore 118511 - 
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S.3.11. Thorium 

Thorium (Th) deposits in the Earth's crust around the world are expected to be three times larg er than those of U, as determined 
from isotope lifetimes and the composition of the accreted material which formed the planet [85-88H However, known resources 
of Th are much smaller and less detailed than those of U, a situation which is the result of the relatively small interest that has been 
given to Th in comparison to U. Therefore, it is to be expected that Th resources increase in size significantly if at some time in 
the future interest grows in Th based nuclear reactors. Although the Th nuclear fuel cycle has been demonstrated several decades 
ago, it has not been used commercially as it involves more safety hazards related to radiation than the U fuel cycle ll8~5ll86ll . The Th 
nuclear fuel cycle is more efficient than that of U and therefore, involves less mass of Th per unit of electricity produced. For similar 
mining costs, Th resources are less expensive per unit of energy, however, the processing of Th into 233 U for fuel preparation has 
not been performed at an industrial scale. Therefore, the cost variable in the Th cost-supply curve is highly uncertain. 

Data for Th resources were obtained from the IAEA y26|. These are provided with much less detail than for U, with four 
uncertainty ranges but only one cost category. Consequently, the strategy of curve-fitting cannot be applied here, and one distribution 
of the hierarchical type per uncertainty category per world region was parametrised in the same way as for fossil resources. The 
cost axis was transformed from a cost per unit of mass to a cost per unit of energy using the efficiency of the Indian experimental 
model reported by Sinha and Kakodkar of 2100 TJ/t 118611 . The curve for the lower boundary of the uncertainty range includes RAR 
only. The most probable curve includes RAR and inferred reserves. The curve for the upper boundary of the uncertainty range 
includes all four categories. 
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S.4. Data tables and figures 
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Figure S.4.1: Curve fits using non-linear least-squares of the two types of cumulative distribution with data from various studies of renewable energy potentials 
previously reported, calculated using the model IMAGE (reproduced from LU d IH). The goodness of these fits are a good indication for which type of 
distribution represents best each type of resource. It can observed that data for abandoned agricultural land is well described by the cumulative distribution for 
nearly identical resources (top left), while the data for rest land is described by the cumulative distribution of the hierarchical type (top right). Data for wind energy 
is well described by hierarchical resource distribution (bottom left), while the data for solar energy is well described by a distribution for nearly identical resources 
(bottom right). 
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■ 1- USA 


HH 


^^■2- Canada 


12 


■ 3- EU-15 


^■13 


1 4- Rest-EU 


^■14 


^^■5- Russia 




^^■6- China 




^^■7- Japan 




^^■8- India 




■ 9- Rest Asia 




^^■10- Oceania 






Region 



Member countries 



USA 

Canada 

EU-15 

Rest Europe 

Russia 
China 
Japan 
India 
Rest Asia 



Oceania 
Brazil 

Rest America 
Africa 
Middle East 



USA 
Canada 

Austria, Belgium, Danemark, Finland, France, Germany, Greece, Ireland, Italy, Luxemburg, Netherlands, Portugal, Spain, 
Sweden, United Kingdom 

Albania, Belarus, Bosnia-Herzegovina, Bulgaria, Croatia, Cyprus, Czech Rep., Estonia, Hungary, Iceland, Latvia, Lithuania, 

Macedonia, Malta, Moldova, Montenegro, Norway, Poland, Romania, Serbia, Slovakia, Slovenia, Switzerland, Turkey, Ukraine 

Russia 

China 

Japan 

India 

Afghanistan, Armenia, Azerbaijan, Bangladesh, Bhutan, Brunei, Cambodia, Georgia, Hong Kong, Indonesia, Kazakhtan, Korea, 

Kyrghizstan, Laos, Malaysia, Mongolia, Myanmar, Nepal, Pakistan, Philippines, Singapore, Taiwan, Tajikistan, Thailand, 

Turkmenistan, Uzbeksitan, Viet Nam 

Australia, New Zealand, Papua New Guinea, pacific islands 

Brazil 

Mexico, Central America, South America excluding Brazil 
Africa 

Barhain, Iran, Iraq, Israel, Jordan, Kuwait, Lebanon, Oman, Palestine, Qatar, Saudi Arabia, Syria, UAE, Yemen 



Table S.4. 1 : Definition of world regions for this paper with member countries for each. 
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Renewable cost-supply curve parameters 





Hydro 


Wind 


Solar 




Hierarchical 






Hierarchical 






Nearly identical 






Region 


A 


B 


Co 


A 


B 


Co 


A 


B 


Co 


Name 


PJ/y 


$/MWh 


$/MWh 


PJ/y 


$/MWh 


$/MWh 


PJ/y 


$/MWh 


$/MWh 


USA 


5 746 


198.40 


1.40 


75 600 


30.19 


145.53 


262 800 


350.03 


620.44 


Canada 


4217 


38.82 


78.06 


43 505 


10.17 


149.28 


51 571 


1542.58 


839.26 


EU-15 


3 066 


22.01 


88.25 


12 009 


65.77 


113.60 


61471 


817.28 


616.21 


R. Eur. 


4 283 


54.13 


38.85 


32 183 


19.38 


142.14 


57 355 


593.61 


718.00 


Russia 


6 595 


82.84 


75.89 


38 774 


65.10 


129.73 


384 448 


428.97 


611.15 


China 


11 502 


26.37 


35.65 


6 057 


175.01 


161.04 


175 509 


241.71 


695.79 


Japan 


846 


25.24 


82.31 


360 


109.24 


161.04 


3 600 


185.69 


840.10 


India 


2 788 


103.23 


0.00 


2 018 


109.62 


208.48 


120 525 


140.01 


483.33 


R. Asia 


7 227 


160.05 


0.00 


18 125 


83.64 


123.77 


282 441 


369.67 


498.17 


Oceania 


779 


71.64 


0.00 


50 410 


59.23 


138.60 


430 421 


141.75 


521.29 


Brazil 


4 738 


18.45 


12.00 


13 248 


23.33 


127.38 


113 386 


310.45 


527.52 


R. Amer. 


7 414 


127.92 


0.00 


22 752 


30.14 


132.96 


160 214 


268.27 


527.44 


Africa 


5 767 


69.16 


64.42 


23 106 


143.77 


161.04 


974 259 


157.16 


456.31 


Mid. East 


1 094 


304.10 


0.00 


7 200 


109.62 


208.48 


306 000 


172.80 


437.87 


Total 


66 061 


345 348 


3 384 000 



Table S.4.2: Table of cost-supply curve parameters for each region for hydro, wind and solar power. 



Renewable cost-supply curve parameters 





Wave 


Tidal 




Hierarchical 






Hierarchical 






Region 


A 


B 


Co 


A 


B 


Co 


Name 


PJ/y 


$/MWh 


$/MWh 


PJ/y 


$/MWh 


$/MWh 


USA 


496 


32.46 


199.44 


145 


89.18 


303.33 


Canada 


5030 


32.46 


199.44 


757 


89.18 


303.33 


EU-15 


1525 


32.46 


199.44 


287 


89.18 


303.33 


R. Eur. 


2442 


32.46 


199.44 


333 


89.18 


303.33 


Russia 


937 


32.46 


199.44 


743 


89.18 


303.33 


China 


361 


32.46 


199.44 


49 


89.18 


303.33 


Japan 


741 


32.46 


199.44 


101 


89.18 


303.33 


India 


174 


32.46 


199.44 


89 


89.18 


303.33 


R. Asia 


2843 


32.46 


199.44 


396 


89.18 


303.33 


Oceania 


1536 


32.46 


199.44 


238 


89.18 


303.33 


Brazil 


186 


32.46 


199.44 


25 


89.18 


303.33 


R. Amer. 


1303 


32.46 


199.44 


253 


89.18 


303.33 


Africa 


1027 


32.46 


199.44 


140 


89.18 


303.33 


Mid. East 


309 


32.46 


199.44 


42 


89.18 


303.33 


Total 


18 910 


3600 



Table S.4.3: Table of cost-supply curve parameters for each region for wave and tidal energy. 
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Primary biomass cost-supply curve parameters 





Al 


A2 


— : 

Keg ion 


A 

A 


r» 
D 


Co 


Type 


A 

A 


r> 
D 


Co 


Type 


Name 


PJ/y 


i/MWh 


VJVlWn 




PJ/y 


i/MWh 


& f\ A~\\ 71^ 

VJVlWn 




USA 


53 082 


3.93 


3.60 


2 


33 082 


6.42 


3.10 


2 


Canada 


18 000 


1.64 


3.60 


2 


12 000 


1.50 


3.60 


2 


CTT K 
DU-1J 


11 305 


1.62 


6.15 


2 


11 305 


0.93 


6.18 


2 


Rest Europe 


11 855 


0.43 


6.28 


2 


10 781 


1.45 


3.44 


2 


Russia 


126 017 


1.73 


3.60 


2 


67474 


1.96 


3.16 


2 


China 


107 100 


11.26 


5.04 


1 


23 322 


9.92 


7.20 


2 


Japan 


115 


0.69 


7.20 


2 


44 


9.98 


14.40 


2 


India 


25 890 


3.37 


3.00 


2 


13 756 


2.32 


2.86 


2 


Rest Asia 


12 921 


0.91 


3.48 


2 


8 605 


0.35 


3.52 


2 


Oceania 


55 625 


2.01 


3.14 


2 


34 477 


3.31 


2.52 


2 


Brazil 


77 751 


1.33 


6.54 


2 


22 070 


5.12 


3.81 


2 


Rest America 


27 629 


2.06 


5.45 


2 


7 310 


3.58 


3.44 


2 


Africa 


139 246 


4.32 


1.68 


2 


53 240 


6.78 


-0.86 


2 


Middle East 


13 011 


25.38 


6.96 


1 


8 011 


13.08 


7.20 


2 


Total 


679 548 


305 477 




Bl 


B2 


— : 

Region 


A 


B 


Co 


Type 


A 


B 


Co 


Type 


Name 


PJ/y 


$/MWh 


$/MWh 




PJ/y 


$/MWh 


$/MWh 




USA 


36 082 


0.85 


5.00 


2 


49 082 


2.07 


3.60 


2 


Canada 


14 000 


0.84 


3.60 


2 


13 000 


0.94 


3.60 


2 


T^T T 1 ^ 


7 268 


3.29 


5.93 


1 


12 921 


4.77 


4.48 


1 


Rest Europe 


9 844 


2.42 


3.59 


1 


12 178 


0.45 


5.70 


2 


Russia 


87 319 


0.99 


3.60 


2 


77 396 


0.94 


3.60 


2 


China 


77 179 


1.57 


3.60 


2 


46 261 


5.51 


7.20 


2 


Japan 


115 


0.69 


7.20 


2 


225 


0.69 


7.20 


2 


India 


13 756 


3.35 


2.92 


2 


6 289 


5.24 


3.60 


2 


Rest Asia 


5 075 


0.72 


3.49 


2 


5 349 


1.80 


3.85 


2 


Oceania 


35 279 


0.96 


2.81 


2 


30 329 


0.87 


3.09 


2 


Brazil 


56 539 


4.76 


4.44 


1 


38 863 


2.29 


6.03 


2 


Rest America 


18 841 


3.86 


3.00 


2 


10 517 


12.20 


3.60 


1 


Africa 


81 240 


3.36 


2.50 


2 


15 237 


4.47 


1.17 


2 


Middle East 


4 011 


4.99 


7.20 


2 


3011 


5.49 


7.20 


2 


Total 


446 548 


320 658 



Table S.4.4: Table of cost-supply curve parameters for biomass primary energy resources for four SRES scenarios Al, A2, Bl and B2. 



Geothermal energy cost-supply curve parameters 





Direct Use 


Electricity 




In belt 






Out of belt 




In belt 






Out of belt 








Hierarchical 






Nearly identical 




Hierarchical 




Nearly identical 




Region 


A 


B 


Co 


A 


B 


Co 


A 


B 


Co 


A 


B 


Co 


Name 


PJ/y $/MWh 


$/MWh 


PJ/y 


$/MWh 


$/MWh 


PJ/y 


$/MWh 


$/MWh 


PJ/y 


$/MWh 


$/MWh 


USA 


74 


7.22 


94.34 


138 


81.44 


118.44 


1 361 


20.62 


144.96 


2 528 


63.62 


255.71 


Canada 


9 


7.22 


94.34 


83 


81.44 


118.44 


253 


20.62 


144.96 


2 280 


63.62 


255.71 


EU-15 


1 


7.22 


94.34 


89 


81.44 


118.44 


14 


20.62 


144.96 


1458 


63.62 


255.71 


R. Eur. 


12 


7.22 


94.34 


27 


81.44 


118.44 


123 


20.62 


144.96 


524 


63.62 


255.71 


Russia 


9 


7.22 


94.34 


177 


81.44 


118.44 


275 


20.62 


144.96 


5 226 


63.62 


255.71 


China 


44 


7.22 


94.34 


103 


81.44 


118.44 


809 


20.62 


144.96 


1 887 


63.62 


255.71 


Japan 


15 


7.22 


94.34 





81.44 


118.44 


144 


20.62 


144.96 





63.62 


255.71 


India 


2 


7.22 


94.34 


37 


81.44 


118.44 


39 


20.62 


144.96 


743 


63.62 


255.71 


R. Asia 


115 


7.22 


94.34 


90 


81.44 


118.44 


1 498 


20.62 


144.96 


1482 


63.62 


255.71 


Oceania 


20 


7.22 


94.34 


80 


81.44 


118.44 


222 


20.62 


144.96 


1 647 


63.62 


255.71 


Brazil 


5 


7.22 


94.34 


96 


81.44 


118.44 


102 


20.62 


144.96 


1 942 


63.62 


255.71 


R. Amer. 


143 


7.22 


94.34 


91 


81.44 


118.44 


2 030 


20.62 


144.96 


1 552 


63.62 


255.71 


Africa 


54 


7.22 


94.34 


210 


81.44 


118.44 


761 


20.62 


144.96 


4 179 


63.62 


255.71 


Mid East. 


11 


7.22 


94.34 


49 


81.44 


118.44 


188 


20.62 


144.96 


901 


63.62 


255.71 


Total 


514 


1 269 


7 818 


26 349 



Table S.4.5: Table of cost-supply parameters for geothermal energy, for both direct use of heat and electricity production. 
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Oil Mtoe 





Crude Oil 




Oil Shales 


Oil Sands 






Extra Heavy Oil 




Region 


[13J 




[15J 


[151 






[15J 






Name 


Reserves 


Resources 


Resources 


Reserves 


Resources 


Additional 


Reserves 


Resources 


Additional 


USA 


3 863 


10 000 


536 931 





5 429 


2 388 


3 


379 


4 


Canada 


667 


2 400 


2 192 


24 909 


227 189 


355 828 











EU-15 


1 193 


1 545 


13 248 


31 


276 





24 


1 928 





Rest Europe 


1 155 


3 530 


4411 





1 





5 


51 





Russia 


10 436 


16 400 


35 470 


4 147 


39 034 


7 505 


1 


25 





China 


2 018 


2 300 


47 600 





233 





110 


1 168 





Japan 


6 


10 























India 


792 


400 























Rest Asia 


9 218 


10 535 


3 988 


6 203 


55 945 





18 


1 163 





Oceania 


595 


1 100 


4 534 




















Brazil 


2 450 


5 000 


11 734 




















Rest America 


8 996 


9 588 


60 





136 





8 476 


270 637 


27 704 


Africa 


17 277 


15 485 


23 317 


263 


2 364 


6 778 


7 


66 





Middle East 


102 366 


21 170 


5 792 




















Total 


161 031 


99 463 


689 277 


35 552 


330 607 


372 500 


8 644 


275 416 


27 707 



Costs USD2008 / boe [801 





Crude Oil 




Oil Shales 


Oil Sands 






Extra Heavy Oil 






Reserves 


Resources 


Resources 


Reserves 


Resources 


Additional 


Reserves Resources 


Additional 


Upper 


10 


10 


50 


40 


40 


40 


10 40 


10 


Lower 


40 


100 


100 


50 


70 


70 


50 70 


70 



Table S.4.6: Oil resources by world region in units of Mtoe (million tonnes of oil). 

Gas Wn? 





Conv. gas 




Shale gas 




Tight gas 




Coalbed Methane 


Methane Hydrates 


Region 


[13J 




[811 




[13J 




[82J 




[83J 




Name 


Reserves 


Resources 


Reserves 


Resources 


Reserves 


Resources 


Reserves 


Resources 


Reserves 


Resources 


USA 


7 080 


20 000 


17 000 


45 600 


1000 


210 000 


9 700 


2 000 








Canada 


1 754 


7 000 


10 988 


42 198 





7 000 


5 700 


70 800 








EU-15 


2 338 


2 530 


7 024 


28 547 





7 000 


2 802 











R. Eur. 


3 889 


7510 


9 374 


39 734 








1 908 











Russia 


47 578 


105 000 


538 


2 152 





45 000 


17 000 


96 300 








China 


2 455 


10 000 


36 109 


144 463 





9 000 


30 000 


5 100 








Japan 


21 


5 


























India 


1 115 


900 


1 784 


8 213 





1 000 


800 











Rest Asia 


23 951 


22 805 


1 444 


5 834 








1 100 











Oceania 


3 553 


2 450 


11 215 


39 111 





1 000 


8 500 


5 700 








Brazil 


365 


2 000 


7 533 


25 658 





6 000 














R. Amer. 


7 704 


8 858 


47 579 


170 745 




















Africa 


14 753 


16 155 


29 482 


112 206 








800 











Mid. East 


75 358 


35 370 


























International 

















38 000 








300 000 


300 000 


Total 


191 914 


240 583 


180 070 


1 0748 62 


1 000 


324 000 


78 310 


179 900 


300 000 


300 000 



Costs USD2008 / GJ f84] 





Conv. gas 




Shale gas 




Tight gas 




Coalbed Methane 


Methane Hydrates 




Reserves 


Resources 


Reserves 


Resources 


Reserves 


Resources 


Reserves Resources 


Reserves Resources 


Upper 


0.5 


0.5 


3.8 


3.8 


2.6 


2.6 


3.8 3.8 


4.4 4.4 


Lower 


5.7 


5.7 


8.6 


8.6 


7.6 


7.6 


7.6 7.6 


8.6 8.6 



Table S.4.7: Natural gas resources by world region in units of Gm 3 (billion cube meters). 
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Coal Mt 





Hard coal 




Reserves 


Resources 


Region 


Proven 


Probable 


Possible 


Proven 


Probable 


Possible 


Name 




[13] 




[151. 


[13J 




USA 


226 694 








6 691 942 








Canada 


4 346 








187 606 








EU-15 


84 721 








278 420 








Rest Europe 


24 534 


752 


1 862 


254 658 


7 428 


11 422 


Russia 


68 655 








2 730 810 








China 


180 600 








681 600 








Japan 


340 








4 603 


I 988 


7 375 


India 


56 100 








105 820 


123 470 


37 920 


Rest Asia 


54 678 








289 048 








Oceania 


44 627 








1 620 675 








Brazil 


1 547 








6212 








Rest America 


9 960 


4572 


4 237 


20 496 








Africa 


32 546 








58 150 








Middle East 


1 203 








41 203 








Total 


790 551 


5 324 


6 099 


12 971 243 


132 886 


56 717 




Soft coal 




Reserves 


Resources 


Region 


Proven 


Probable 


Possible 


Proven 


Probable 


Possible 


Name 


[15J 


[13] 




[15] 


[13J 




USA 


30 851 








1 398 669 








Canada 


3 108 








17 371 


40 055 


108 995 


EU-15 


44 214 








89 158 








Rest Europe 


40 456 


1 996 


3 124 


275 185 


14 961 


11 581 


Russia 


91 350 








1 371 030 








China 


52 300 








318 000 








Japan 


10 








160 


1 132 


4 074 


India 


4 895 








38 647 








Rest Asia 


30 762 


7 086 


34 070 


387 263 


11 871 


57 198 


Oceania 


37 738 


62 840 


101 100 


46 973 


73 102 


112 300 


Brazil 


4 559 


7 559 


4 575 


6 513 


10 799 


6 535 


Rest America 


5 633 


527 


790 


7 524 








Africa 


180 








338 








Middle East 




















Total 


346 056 


80 008 


143 659 


3 956 831 


151 920 


300 683 




Costs 


USD2008 / 1 


[80] 












Reserves 


Resources 




Proven 


Probable 


Possible 


Proven 


Probable 


Possible 




20 


20 


20 


20 


20 


20 




50 


50 


50 


100 


100 


100 



Table S.4.8: Coal resources by world region in units of Mt (million tonnes of coal). Hard coal includes anthracite and bituminous coal, while soft coal includes 
sub-bituminous coal and lignite. Since there is no clear demarcation between ranks of coal, the limit is put onto the calorific content, and thus coal resources with a 
calorific content higher than 16 500 kj/t belong to the hard coal category (as defined in [ 13]), while coal resources with a lower calorific content belong to soft coal. 
Anthracite can have calorific contents of up to 35 000 kJ/t while lignite can have calorific values as low as 1 1 000 kj/t. 
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Uranium 


t 


[26] 














Region 


Reasonably Assured Reserves (RAR) 


Inferred 


Name 


<40$/kg 


<80$/kg 


<130$/kg 


<260$/kg 


<40$/kg 


<8()$/kg 


<130$/kg 


<260$/kg 


USA 





39 000 


207 400 


472 100 





19 500 


103 700 


236 050 


Canada 


267 100 


336 800 


361 100 


387 400 


99 700 


1 10 600 


124 200 


157 200 


EU-15 





7 000 


20 800 


33 800 








13 500 


110 400 


Rest Europe 


2 500 


39 100 


88 500 


160 000 


3 200 


15 000 


43 850 


109 850 


Russia 





100 400 


181 400 


181 400 





57 700 


298 900 


384 900 


China 


52 000 


100 900 


115 900 


115 900 


15 400 


49 100 


55 500 


55 500 


Japan 








6 600 


6 600 








3 300 


3 300 


India 








55 200 


55 200 








23 900 


24 900 


Rest Asia 


14 600 


326 600 


454 500 


533 500 


29 800 


276 800 


366 000 


474 900 


Oceania 





1 163 000 


1 176 000 


1 179 000 





449 000 


497 000 


500 000 


Brazil 


139 900 


157 700 


157 700 


157 700 





73 600 


121 000 


121 000 


Rest America 





7 000 


11 700 


13 800 





4 400 


10 100 


11 300 


Africa 


93 800 


194 600 


644 100 


663 400 


78 500 


121 800 


260 000 


286 300 


iviiQcue cast 





44 000 


44 000 


44 700 





67 800 


67 800 


69 200 


Total 


569 900 


2 516 100 


3 524 900 


4 004 500 


226 600 


1 245 300 


1 989 750 


2 544 800 


Region 


Prognosticated 


Speculative 


Name 




<80$/kg 


<130$/kg 


<260$/kg 




<130$/kg 


<260$/kg 


Unassigned 


USA 




819 500 


1 169 300 


1 036 950 




858 000 


858 000 


482 000 


Canada 




50 000 


150 000 


150 000 




700 000 


700 000 





EU-15 




7 000 


7 600 


7 600 




50 100 


50 100 


94 000 


Rest Europe 




200 


41 650 


66 650 




6 650 


126 650 


314 000 


Russia 







182 000 


182 000 










633 000 


China 




3600 


3 600 


3 600 




4 100 


4 100 





Japan 







3 300 


3 300 




3 300 


3 300 





India 










63 600 










17 000 


Rest Asia 




377 900 


591 400 


592 900 




1 776 600 


1 806 100 


264 700 


Oceania 




300 000 


300 000 


300 000 










500 000 


Brazil 




73 600 


121 000 


121 000 




121 000 


121 000 





Rest America 




6 600 


23 500 


23 500 




236 700 


236 700 


176 200 


Africa 




49 400 


156 900 


156 900 




25 000 


25 500 


1 112 900 


Middle East 




67 800 


89 000 


89 000 




84 800 


98 800 





Total 




1 755 600 


2 839 250 


2 797 000 




3 866 250 


4 029 750 


3 593 800 


Table S.4.9: Uranium resources (in natural concentration) by world region in units of tonnes. 


Thorium 




t 




[26J 










Region 




RAR 




Inferred 




Indentified 




Prognosticated 


Name 




< 80 USD/kg 




< 80 USD/kg 




< 80 USD/kg 




N/A 


USA 




122 000 




278 000 




400 000 




274 000 


Canada 









44 000 




44 000 




128 000 


EU-15 






















Rest Europe 




54 000 




213 000 




186 000 




164 000 


Russia 




75 000 




112 500 




75 000 







China 






















Japan 






















India 




319 000 




478 500 




319 000 







Rest Asia 






















Oceania 




46 000 




406 000 




452 000 







Brazil 




172 000 




130 000 




302 000 




330 000 


Rest America 









300 000 




300 000 







Africa 




18 000 




127 000 




118 000 




410 000 


Middle East 






















Unassigned 




23 000 




10 000 




33 000 




81 000 


Total 




829 000 




2 099 000 




2 229 000 




1 387 000 




Costs 




USD2008 / kg 


















RAR 




Inferred 




Indentified 




Prognosticated 


Lower 




40 




40 




40 




80 


Upper 




80 




80 




80 




260 



Table S.4. 10: Thorium resources (in natural concentration) by world region in units of tonnes. 
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